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the present study deals with the subject of optimal 
balancing of planar mechanisms when the contribution due to 
vibratory inertia forces of mechanisms with elastic links 
towards’ the shaking force and the shaking moment transmitted 
to the supports is taken into consideration along with the 
contribution due to the rigid body inertia forces. A 
reduction in the shaking force, the shaking moment and the 
input end torque fluctuations is obtained by an appropriate 
reduction in the total mechanism mass, i'his is achieved by 
allowing the cross sections of the links to vary and putting 
constraints on the stresses in the links and on the deflec- 
tions of the strategic points on the mechanism.. A procedure 
to obtain an optimum, design of planar mechanisms is presented 
using different methods of achieving the balance in these 
mechanisms. For a specified aim, which the mechanism has 
to serve, the kinematic synthesis of the mechanism is 



chazLged to arrive at an optimum design from the point of 
view of balance in. the mechanism incorporating the afore- 
mentioned aspects. The effect of redistributing the link 
masses in order to reduce the shaking force and/or the 
shaking moment on the stresses in the links and on the 
contribution of vibratory inertia forces towards the shaking 
force and the shaking moment has been studied and an optimum 
design is obtained taking these effects into consideration. 

The kinetoelastodynamic analysts of the mechanism is 
made to determine the stresses in the links and the vibratory 
inertia forces using the finite element approach of struc- 
tural dynamics. Taking the combined effect of the rigid 
body inertia forces and the vibratory inertia forces the 
maximum shaking force and the maximum shaking moment during 
a complete motion cycle of the mechanism are determined 
using the dynamic force analysis. 

The constrained minimization problem is first converted 
into an unconstrained problem using interior penalty function 
method , Davidson-flet cher-Powell method is then used to 
solve it. Since the evaluation of the stresses in the links 
and the vibratory inertia forces at each design point during 
the optimization process is very time consuming, two approxi- 
mate procedures have been adopted, firstly, during one 
dimensional minimization the gradients of the maximum 



link stresses, maximum shaking force and the maximum shaking 
moment are used to predict them at the neighboring design 
points. Secondly, during the calculation of the abo-^T-e men- 
tioned gradients and the gradient of the penalty function by 
backward finite difference method, the kinetoelast odynami c 
analysis, to obtain the maximum stresses in the links, ttie 
maximum shaking force and the maximum shaking moment, is made 
for the com.plete cycle of mechanism motion only at the base 
point. When each design variable is perturbed- by a sm.all 
aro.ount the kinetoelastodynamic analysis is, however, made 
only upto a certain small interval beyond the position at 
which either the maximum stress in the links or the maximum 
shaking force or the maximum shaking moment (whichever 
corresponds to maximum crank angle) has occurred at the 
base point. Ihese approximations have been obsenjed to be 
reasonably good and have reduced the computational time 
considerably. 

The results obtained show that to achieve a balance 
in mechanism with elastic links in a real sense the contri- 
bution of vibratory inertia forces towards the shaking force 
and the shaking moment must be taken into consideration, 
furthermore reducing the areas of the cross-sections of the . 
links by putting constraints on the stresses in the links 
is very effective in minimizing the shaking force and the 
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shaking moment as wall as the fluctuations in the input enc? 
torque. Thus, the present work represents a unified approach 
to the design of optimally balanced mechanisms with elastic 
links . 



CHAP'iER I 


INTRODUCTION 

Every machine is composed of one or more mechanisms 
performing different tasks. The successful operation of a 
machine depends on the extent the expected operating 
conditions are realized during the process of synthesis 
through various assumptions made at that stage. Kinematic 
criteria were the major considerations in the design of 
mechanisms in the long past. With the advent of machines 
moving at higher speeds the consideration of .dynamic 
criteria for synthesis of mechanisms have become essential. 
The dynamic criteria for synthesis such as fluctuations 
in input speed , balance of mechanism, and stresses due to 
dynamic forces were based on ’che assumption of rigid links 
till recent past. However, during the last decade, research 
has shown that the results predicted with the assumption of 
rigid links are far from accurate due to the elasticity of' 
the links,' specially at high speeds or when the links are 
massive. Thus, there is a need to put the design techniques 
on more rigorous basis by developing new synthesis proce- 
dures that take into account the elastic properties of the 
mechanism. The present study is an attempt in this direc- 
tion. 



When a mechanism moves, occupying various positions in 
space, the velocity and acceleration of the links vary during 
a complete cycle of motion of the mechanism. thus, there are 
varying forces and moments due to rigid body inertia of the 
links transmitted to the machine foundation which results in 
undesirable vibrations of the machine as well as its supports 
As the speed of the mechanism increases these unbalance 
■forces and moments transmitted to the foundation pose a much 
greater problem. Thus, it becomes necessary to reduce these 
forces and moments. Several methods are available for this 
purpose. However, in all these methods the consideration 
is limited to mechanisms with rigid links. In actual prac- 
tice, specially at higher speeds, the links of a mechanism, 
due to their elasticity, vibrate about some mean position 
under the action of rigid body forces and external forces 
present on the mechanism. These vibrations give rise to an 
additional acceleration field besides the acceleration field 
produced by the gross rigid body motion of the links. The 
acceleration field resulting from the vibration of the links 
develops additional inertia forces which may be called as 
kineto-elastodynami c inertia forces (referred to as KBD 
inertia forces hereafter). The order of magnitude of OD 
inertia forces depends upon the mass distribution of the 
links as well as the magniuude of the acceleration field 
resulting from the vibration of the links. If any of these 
two is large the resulting ESD inertia forces will also be 
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high. As such, the varying forces and moments transmitted 
to the machine foundation due to the FRT ) inertia forces -will 
also be substantial. Therefore, while considering the 
problem of reducing or balancing the unbalance forces and/or 
moments transmitted, to the machine foundation due to rigid 
body inertia, one should also take into account the effect 
of KBD inertia. furthermore, while attempting to reduce or 
balance the forces and moments transmitted to the foundations 
due to rigid body inertia, either by redistribution of the 
link masses or by changing kinematic synthesis, no effort 
has been made to observe the effect of these changes on the 
stress distribution in the links. Also no attempt has been 
made to reduce the unbalance forces and moments by a reduc- 
tion in the total mechanism mass by reducing the areas of 
cross-sections of the links, which is not possible without 
the inclusion of the analysis of stresses in the links. 

Thus, for the rational dynamical design of mechanisms 
with regard to minimization of the shaking force and/or 
the shaking moment one should also include the effect of 
KED inertia and keep the stresses in the links and the 
deflections of the vital points on the mechanism within 
prescribed limits. Moreover, efforts should he made to 
develop approximate methods to obtain the solution of such 
problems with reduced computer time without any significant 
loss of accuracy. 
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A brief survey of the publicatins relevant to this work 
is presented below. 

1.1 A Survey of Past inForks 

1.1.1 Mechanism Balancing 

'I'he subject of djmamics of mechanisms deals with the 
motion of mechanisms under the action of actuating forces and 
torques and also the forces and the torques produced by the 
motion of the mechanism. Mechanism balancing involves the 
design or modification of a mechanism to control the force 
and the torque levels produced by its motion. Ihe mechanism 
balancing falls under two major categories, namely, (i) input 
torque balancing and (ii) shaking force and shaking moment 
balancing. 

1.1.1(a) Input I'orque Balancing 

Fluctuations in input speed and torque in a mechanism 
occur since the ehergy requirements vary during a complete 
cycle of motion of the mechanism. I'o reduce or minimize these 
fluctuations many researchers have put forx-fard various 
procedures. Some important ones are briefly discussed below. 

Sherwood [1] attempted to make the total kinetic 
energy of a mechanism constant (or approximately constant:) 
during its entire cycle of motion. I'his resulted in constant 
input shaft Speed or reduced, variation in the input shafb 



5 


speed. Thus, for no load or constant output load the input 
end torque fluctuations vjere minimised. This ’iso-energetic' 
condition was achieved by the optimum distribution of link 
masses in the mechanism. For a four bar linkage a certain 
distribution of coupler mass made the total kinetic energy 
approximately constant. For a drag link mechanism the 
results were very good, whereas for a crank- rocker mechanism, 
the variation in the total kinetic energy was reduced to 
half only, 

Sherwood and Hockey [2] presented a general method of 
optimization of the mass distribution in the links of a 
mechanism using dynamically similar system which minimised 
the fluctuation in the kinetic energy of the mechanism run- 
ning at the constant angular input speed. Using the properfy 
that any link in plane motion might be represented by four 
variable point masses at four almost arbitrarily chosen 
points fixed in the link, the links of the mechanism were 
represented hy the above mentioned dynamically equivalent 
set of point masses and an optimum mass distribution was 
worked out giving minimum fluctuation in the kinetic energy. 

Skreiner [3J studied the fluctuation of the input 
torque and demonstrated that it could be red.uced by the 
addition of a spring in the system of linkages. Assuming 
the links to be rigid he performed the dynamic motion 
analysis for the resulting one degree of freedom system- 
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Applying ihe principle of virtual work the equation of motion 
was established. He determined the variation in input speed 
and the pin forces in one case by varying the inertia of the 
input shaft and in another case by using a spring with dif- 
ferent attachm.ent points to the frame. He showed that ft 
was possible to reduce fluctuation in the input speed without 
much increase in the pin force by the use of a suitable 
spring . 

Hockey [4j presented an improved technique for obtain- 
ing the optimum mass distribution of the links to reduce the 
fluctuation of kinetic energy in planar mechanisms over their 
previous procedure [2] so as to place the concept on a more 
rigorous mathematical footing. He concluded that the fluc- 
tuation of kinetic energy of a mechanism could be efficiently 
reduced, by proper distribution of link masses and it was 
claimed that his new method of determining this 'proper' 
distribution was superior to the previous method. He 
reported, that this method of mass distribution was more 
efficient than the use of a flywheel at the input shaft . 
However, he pointed out that though this technique would 
reduce the fluctuations of kinetic energy and input speed 
but other dynamic effects such as force and moment unbalance 
might worsen. 

Hockey [5] ’took up the problem of minimizing fluctua- 
tions of input torque resulting not only due to the inertia 
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of the moving links bux also due to fluctuating external 
loads. The previous works, which were aimed at reducing 
the fluctuation of kinetic energy in a mechanism in the 
absence of external loads and friction, also resulted in. 
the reduction in the fluctuation of torque required at xhe 
input shaft to maintain constant speed.- He developed a 
method of dynamic synthesis using system of variable point 
masses as outlined in his earlier work [4] to include 
mechanisms which had resisting loads on their links. The 
mass in the links of the mechanism was so distributed that 
the fluctuation of the torque required to overcome the 
external loads was counteracted by the fluctuation of torque 
required to overcome the effects of inertia. It was shown 
that the mechanism might be completely force balanced by 
using a method due to Berkof and Lowen [6] vjithout affect- 
ing in any way the minimized fluctuation of input shaft 
torque. 

1.1.1(b) Shaking force and Shaking Moment Balancing 

Lowen and Berkof [7] presented a detailed survey of 
the force and moment balancing literature for mechanisms 
with rigid link of constant mass. Eleven translations from 
the German and Russian literature in English were also 
included. Therefore in this survey, only references of 
relatively recent origin have been discussed . 



s 


Han [8] developeci an analytical method for the optimal 
balancing of the shaking force and shaking moment of forces 
of a planar mechanism, driven by a constant speed shaft. He 
used a single balance weight and determined its mass moment 
of inertia and phase angle which gave minimum disturbing 
action through a complete revolution of the driving shaft. 

He applied Legendre's principle of least squares of the 
compromised expression ..for the deviations of the shaking 
load on the machine supports produced by the resultant 
shaking force and shaking moment to obtain the optimum set 
of values for 'the balance weight. He also suggested an 
iterative procedure for determining the mass moment of 
inertia and phase angle of the balance weight when the 
input angular speed was not constant. However, his method 
was not suitable for application to every machine in general. 
Use of a single mass rotating at the crank speed is helpful 
only when the force vector is predominently the first 
harmonic. for other mechanisms, specially in which the 
center of mass of the parts is far from circles, ellipses 
and straight lines, balance .masses rotating at constant 
speed is of not much use. 

Berkof and Lowen [6] developed the "method of linearly 
independent vectors" to obtain the complete force balance in 
planar linkages . i'he position of the total mechanism center 
of mass was written in such a way that the coefficients of 
the time dependent terma could be equated to zero, thereby 
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making the total center of mass stationary. I’hus the 
resultant of all dynamic forces reduced to zero for all 
positions of the mechanism. 'I'he relationships obtained hy 
equating the co -efficients of the time dependent terms to 
zero were used to determine the necessary distribution, of 
link masses. Using this method the shaking force was 
completely eliminated but it caused the bearing reactions, 
the shaking moment and the input torque variations to 
increase. 

Later the same authors [9] presented a least squares 
optimization theory which minimised the shaking moment of a 
fully force balanced inline four bar linkage. Using the 
angular momentum principle an expression was derived for the 
shaking moment of a force balanced four bar linkage. Next, 
a periodic function of the same general type as the shaking 
moment was formulated. I'his periodic function was made to 
deviate least from zero in root -mean-square sense, to yield 
the optimum relationship between the coefficients of the 
periodic function. fhe coefficients of the moment expression 
depend only on the link niass distribution. l*he optimum 
coefficient ratios are functions of link lengths only. 

lowen and Berkof [10] demonstrated how their above- 
mentioned theory could be applied to obtain a linkage with 
an optimum shaking moment characteristic. According to 
these authors, the theory was applicable to most linkages. 
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but for practical reasons the field of application was 
restricoed to standard linkage configurations which were 
generally used in practice. A matching procedure was 
developed to determine all these families of mechanisms for 
which the standard configurations represent the optimum 
mass distribution. I'his method required the construction 
of moment optimization graphs which were used to match the 
theoretical optimum coefficient ratios of the moment equa- 
tion with the actual attainable ratios. 

Kaufman and Sandor [11] extended the 'method of 
linearly independent vectors' [6] to spatial linkages and 
obtained the complete force balance of spatial linkages such 
as the RSSR and RSSP mechanisms . 

iepper and Lowen [12] gave a contour theorem that 
differentiated between mechanisms which could be fully force 
balanced and those which could not. In addition, the method 
of linearly independent vectors [6] has been generalized 
and it was shown that the 'apparent' minimum number of 
counterweights producing full force balancing by internal 
mass redistribution alone equals n/2 where n was the number 
of links in the mechanism. 

Berkof [13] obtained the complete shaking moment; 
balance of force-balanced inline four bar linkage using 
two circular masses moving in opposite directions to that 
of the crank and the follower, respectively* with the help of 
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gears. The moment of inertia of the masses were so chosen 
that the inertial torques produced by them balanced, the con- 
tribution due to the crank and the follower assembly towards 
the shaking momenc respectively. The contribution of the 
coupler Ixnm towards the shaking moment was made zero by 
making the coupler a physical pendulum, This procedure, 
however, resulted in increased bearing reactions and increase 
fluctuations in the input end torque. 

Tepper and lowen [14J presented a method, wherein the 
ms shaking force of a general four bar linkage running at 
constant inpuc. shaft speed was optimized while the IMS ground 
bearing reaction were limited to lie between those corres- : 
ponding to the unbalanced mechanism and those of the fully 
force-balanced mechanism. The optimization was done using a ^ 
single counterweight as well as two counterweights. The 
authors claimed that the two -counterweights technique was 
capable of producing considerably lower bearing forces than i 
the single-counterweight technique for comparable values of i 
the EMS shaking force. However, the authors did not make 
any attempt to limit either the EMS input moment or the EMS i 
shaking moment. The effect on coupler bearing forces was 
also not considered. ’ ^ 

Conte, George, Mayne and Sadler [15] used nonlinear 
programming techniques to combine the kinematic ' synthesis 
and dyaamic design of four-bar mechanisms. The arbitrary 
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parameters of kinematic synthesis were determined to satisfy- 
the performance criteria related to dynamic force, as well as 
kinematics. In their examples, the authors obtained optimum 
crank -rocker path generating mechanisms with prescribed 
timing, using various objective functions such as shaking 
force, shaking moment, input torque fluctuation and bearing 
reactions for minimization. The optimization procedure 
resulted in a significant improvement in these dynamic 
characteristics. Minimizing the total mechanism mass, the 
authors demonstrated that instead of an improvement, all the ; 
dynamic performance properties became worse compared to that : 
of the original linkage, thereby indicating that minimization I 

of the total mechanism mass was not fruitful. ; 

/ 

Wiederrich and Roth [16] deternined condition for 
reducing the angular momentum fluctuations transmitted to | 
the frame of a completely force-balanced four bar linkage. 

The authors claimed that their approach led to simple design ] 
equations for determining the inertial properties of the 
links for good momentum balance, as compared to the computer : 
based search techniques required in earlier methods. 

Elliott and Tesar [17] extended the "method of linearly 
independent vectors" of- Berko f and Lowen [6] to shaking 
moment and driving torque functions. They defined four mass : 
parameters for each link of. the mechanism. These were the | 
mass, the radius of gyration (or mass moment of inertia) and } 



the centroid coordinates. Thus, they cons'idered twelve mass 
parameters in a four bar linkage and by adjusting as many of 
these parameters as possible they obtained a balanced device, 
for this purpose they described the shaking moment as the 
time rate of change of total momentum about the input crank 
pivot and finally put it into a form similar to that of the 
shaking force. By equating the coefficient of time dependent 
terms to zero a relationship between various parameters was 
obtained. But it required two rotary inertias to te nega- 
tive, which was achieved by adding counter rotating masses 
to the system. Thus using their general approach they 
completely eliminated the shaking force and the shaking 
moment with the addition of physical negative mass for all 
the four-link mechanisms; the four bar, the slider crank, 
the inverted slider crank and the oscillating block mecha- 
nisms. Prom a similar formulation for driving torque it 
was possible to meet exactly specified values of driving 
torque upto three positions or in the least -square sen8„- 
for the full cycle while maintaining shaking force balance. 

Walker and Oldham [18] presented a procedure of 
force balancing multidegree of freedom, mult ibar linkages . 

They derived three equations from which the force balance 
conditions were obtained. The linkage could have both 
revolute or prismatic joinus. However, there were force 
balance restrictions with prismatic joints. 
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Carson and Stephens [19] extended the theory developed 
by Berkof and Lowen [9,10] to problem of .force and EJViS moment 
balancing by redistribution of masses when all links length 
ratios were specified. fhe criteria for distributing the 
link mass involved selecting numerical values for eight 
unknown mass parameters to satisfy three equations. Selec- 
tion of five parameters arbitrarily and calculating the 
remaining three might result in physical impossibilities 
such as imaginary radius of gyration. To overcame this 
difficulty the authors developed feasible design spaces and 
monographs showing feasible parameter design spaces and 
parameter interrelationship. 

Using the ’method of linearly independent vectors' 
Balasubramanian and Bagci [20] developed design equations 
for complete shaking force balancing of commonly used planar 
Stephenson’s and Watt’s type 6 R 6-bar and 6-bar slider- 
crank regular force transmission mechanisms. It was aimed 
to provide counter-balancing masses only on the shafts 
supported by the fixed frame whenever possible to avoid 
the worsening of the unbalanced shaking moment. 

Walker and Oldham [21] further extended their earliei? 
work [18] and presented a method which indicated whether a 
linkage could be fully force-balanced using counterweights 
alone. A formula was given, to determine the minimum, number 
of counterweights needed for a full force-balance- A criterio 



was presented to select an optimum counterweight set so 
that the increase in bearing reactions, driving torque and 
out of balance couples were checked. 
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Bagci [22] presented a me chod of complete force- 
balancing of planar mechanisms with force transmission 
irregularities by adding idler loops and using the method 
of linearly independent mass vectors. A mechanism contain- 
ing a link or group of links that had connection to frame 
of the Hiachine through pairs all permitting linear freedom 
was termed as irregular force transmission mechanism. The 
component of shaking force contributed by these links could 
not be balanced by means of distributing the masses of the 
existing links of the mechanism. However by the addition 
of idler loop consisting of a. set of two binary links having 
all revolute pairs to the original mechanism it was made 
possible to completely force balance the mechanism. 

1.1.2 Deflection and Stress Analysis 

Erdman and Sandor [23] presented an excellent review 
of the work carried out oia the deflection and sti’ess analysis 
in the mechanisms due to elastic links till about 1971 • A 
considerable progress has been made in, this direction since 
then by various researchers- Therefore these latter works 
are briefly presented in this section. 



Winfrey [24] used finire element method to determine 
the longitudinal, transverse and torsional displacements of 
a general mechanism. Each link was modelled to have six 
elastic and three rigid body coordinates. Using a connec- 
tivity matrix the element matrices were assembled into 
system matrices. Using conservation of momentum during the 
free vibration of the mechanism a transformation matrix was 
obtained for the elimination of the rigid body degrees of 
freedom. Use of this transformation matrix made the stiff- 
ness matrix non-singular and hence the eigenvalue problem 
in the transformed coordinates could be solved. With the 
help of the modal matrix thus, obtained, and by assuming the 
damping ma.trix to be proportional to the mass matrix, a set 
of uncoupled differential equations in normal coordinates 
was obtained, I’he solution of these equations was deter- 
mined and converted back into the original system coordinate 
from the normal coordinates hy using first the modal matrix 
and then the transformation matrix. I'he final displacement 
and velocity determined at the end of current interval of 
time were used as initial conditions for the next interval. 

Brdman, Sandor and Oakberg [25j used the flexibility 
approach of structural analysis to determine the elastic 
deflection of entire mechanism system. The effect of 
elastic accelerations was also included in determining 
the load vector. Eor this purpose a kineto-elastodynamic 
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stretch rotation operator was defined which was used to 

% 

determine the total acceleration of the mechanisin including 
elastic accelerations. Starting from these accelerations' 
the inertia forces were determined and deflection analysis 
was made. The total accelerations were then again deter- 
mined and compared with the old ones. If they were not 
within acceptable tolerance then above procedure was repeated 
The crank was assumed as a cantilever beam which converted 
the mechanism into structure. 

Erdman, Imam and Sandor [26] extended the above 
analysis and used the method of 'dynamic equivalence system' 
to take into account the elastic inertia forces. It was 
indicated that in certain cases the elastic accelerations 
were comparable with the rigid body accelerations. 

Imam, Sandor and Kramer [27] analysed planar mecha- 
nisms with elastic links using stiffness approach of the 
finite element method . Permutation vector method was used 
to assemble the element nBtrices into system matrices. To 
reduce the computational time the rate of change of eigen- 
values and eigenvectors was employed. The rotation at the 
support end was not included in the system coordinate 
thereby making the crank a cantilever beam and reducing 
the mechanism, to a structure. Thus the true boundary 
conditions of the mechanism were not met during the analysis. 
To determine the stresses, using deflections obtained from 
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the above analysis, a pol 3 momial expansion method based on 
the theory of strength of materials [28] was used. 

Imam and bandor [29] used these methods of kineto- 
elasto dynamic design of elastic mechanisms to minimise the 
mass of the mechanism with constraints on stresses and 
deflections. I'he optimization problem was formulated in 
terms of stepwise linear programming. Simplex method was 
used to .obtain a linkage with minimum mass with linearized 
constraints on deflection and stresses. 'I'o keep the error- 
caused by linearization of constraints within limits the 
range of variation of design variable was restricted, to a 
region about the initial design point by the help of side 
constraints. Later the authors [-30] applied nonlinear 
programming technique to meet the same goal. Both quasi- 
static and rate of change of eigenvalues methods were used 
to determine the elastic deflections. 

Sadler and Sandor [31] presented a lumped mass approach 
to analyse the vibration and stresses in a slider crank 
mechanism. They [32] further extended their method of 
analysis to a planar four-bar linkage with three, elastic 
links. Applying d ' Alem.bert ' s principle to each lumped 
mass the Luler's equation of elastic curve was obtained, 
wMch was then solved using- finite difference method to 
obtain the transverse deflections. The additional relative 
normal and tangential accelerations were taken into account 



and the pin forces were obtained applying force and moment 
balance to each link at their instantaneous positions. The 
additional rigid body motions of the coupler and the follower 
due to elastic deformations of the crank were also considered 
The longitudinal deformations due to axial forces and the 
foreshortening of the links due to bending were neglected. 

Alexander and Lawrence [33,34] were the first to 
obtain experimentally the strain histories of the coupler 
and rocker midpoints of a four bar linkage with elastic 
links. These results were compared with the analytical 
results obtained by using the stiffness method presented by 
Imam, Sandor and Khamer [27] • To represent the continuous 
displacements of the individual links in terms of their 
nodal displacements a two term sine series was assumed. 
Runge-Kutta method was used to solve the linear coupled 
equations. There was reasonable agreement between the 
experimental and analytical results. 

Nath [35] used the basic concept of the finite ele- 
ment technique to analyse elastic mechanisms. He included 
many considerations which were not considered in the earlier 
works. The effect of subdividing the links on the solution 
accuracy was studied. To account for the dynamic behaviour 
of the rigid body inertia forces, the harmonic analysis of 
these forces was made and used for the vibration analysis . 

A novel method of making rigid body analysis by the use of 



a trajisformation matrix was presented, ihe effect of addi- 
tional acceleration terms and the rigid body axial forces 
on the dynamics of the mechanism was observed. A method to 
obtain steady state displacements for all configurations of 
the mechanism by means of a single expression was developed 
for a restricted class of mechanisms incorporating all the 
dynamic factors. However, the validity of this method was 
not well demonstrated. 

Sutherland [36] adopted, a different technique (unlike 
flexibility or stiffness approach of structural dynamics) 
by deriving ohe Euler-Lagrange equation of .motion of a 
purely elastic four bar linkage using the assumed mode 
analysis approach [37]* I'he form of the solution for the 
member deflections was discussed and the governing non- 
dimensional parameters were identified, I'he results were 
determined in a non-dimensional form. These results were 
compared with the results of physical experimentation tc 
establish the validity of the mathematical model. 

G-olebiewski and Sadler [38] compared analytics,! and 
experimental dynamic bending stresses in a elastic connec- 
ting rod of a slider-crank mechanism with a rigid crank. 

The differential equation including viscous damping was 
derived by way of lumped parameter approach using d'Alembert 
principle and the Euler-Bernoulli beam theory, and solved 
numerically. Experimental data were obtained using strain 
gauges . 



Tiiompson. and Barr [30] applied variational procedure 
for setting up the equations describing the elastodynamic 
motion of planar linkages in which all the members were 
considered to be elastic. By permitting independent 
variations of stress, strain, displacement and velocity 
parameters for each link approximate equations of motion, 
boundary and compatibility conditions for the complete 
mechanism were systematically constructed. The procedure 
modelled the physical situation by allowing the simulta- 
neous displacement of all the members from the rigid linked, 
configuration, each displacement being the sum of rigid 
body and deformation components. The authors claimed that 
the following special features could be incorporated with 
this procedure. (i) The links might be made of different 
materials or even anisotropic materials. (ii) Linear or 
nonlinear strain displacement relationships could be taken. 

(iii) A particular link might be considered as rigid 
while the other members of the mechanism were flexible. 

(iv) The effect of viscous torque at each joim; could be 
incorporated, 

Jandarasits and Lowem [40] had analysed the trans- 
verse elastic behaviour of a count erweighted rocker link 
with an endmass. The crank and coupler were taken as rigid. 
The Hamilton's integral and method of Kantorovich were 
used to obtain the linearized decoupled Hill's equations, 
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which were solved to furnish the time portion of the solu- 
tion. 'Ilie space portion of the solution was determined 
using normal modes of free vibration of the complex link. 

An elastic mechanism, constraint equation, relating the 
elastic link angles to the foreshortening of the elastic 
link, provided the . necessary auxiliary conditions. Deter- 
mination of the stability boundaries and steady--sxat e solu- 
tion of the Hill's equations with and without damping was 
presented. In a companion paper [41] they applied their 
technique to an example mechanism and the results were 
compared with those of the experiment conducted by them. 

I'hey confirmed that a good qualitative as well as quantita- 
tive agreement existed between analytical and experimental 
results. 

1.2 Objectives of the Present Work 

Prom the above survey of the past research work, it 
is observed that the efforts aimed at eliminating or reduc- 
ing the shaking force and/or shaking moment transmitted to 
the machine foundations, were based on the consideration of. 
the effect of inertia forces, produced due to accelerations 
resulting from gross rigid body motion of the mechanism, 
towards the unbalance in the mechanism. VJhereas in a real 
situation, specially when the mechanism runs at high speed and 
has elastic links, there are inertia forces arising ouf of 
the accelerations due to 'the vibrations of the links., I'he 



contrilDution of these vibratory inertia forces towards the 
shaking force and the shaking moment, to the best of author's 
knowledge, has not been taken into account in earlier works. 
One may think of making the cross-sections of the links 
heavier in order to reduce their vibrations. But a, careful 
study can reveal that xhis may further worsen the situation 
due to the increased rigid body inertia forces. On the ouher 
hand, a reduction in the cross-sectional areas of the links 
results in a reduction in the rigid body inertia forces 
thereby red.ucing the shaking force and the shaking moment 
transmitted to the supports. 'Ihis can be adopted if the 
stresses in the links are determined at each stage and ke.pt 
within allowable limits. • However, a reduction in the cross- 
sectional areas of the links increases the vibratory accelera 
tions of the links. As such, tho.ugh the inertia of the links 
is decreased, the vibratory inertia forces may be higher. 
Furthermore , the stresses in the links have not been taken 
into consideration while redistributing the link masses to 
achieve a reduction in the unbalance shaking forces or 
shaking moments. In fact, the stress distribution is 
changed when the link masses are redistributed by placing 
counterweights. 

The major objectives of the present work may be 
stated as follows: 



(i) The shaking force and the shaking moment 
resulting from vibratory inertia forces will be taken into 
consideration to design a mechanism with optimal balancing 
characteristics. 

(ii) The method of changing the kinematic design to 
achieve a reduction in the shaking force and the shaking 
moment [15] shall be adopted to design an optimally balanced 
mechanism with constraints on stresses in the links and/or 
constraints on deflection of strategic points on the mechanism. 

(iii) The method of redistributing link masses to 
reduce the shaking forces and the shaking moment will be 
used with constraints on stresses and/or on deflections to 
obtain optimum design of a mechanism. 

(Iv) The areas of cross-sections of the links x^ill 
he varied and its effect on the balancing oharacteristlcs 
will he investigated. 

(v) The determination of the maximum stresses in 
the links and the KBD inertia forces is a time consuming 
process, therefore the investigation of problems involving 
optimization with KEID inertia forces and with stress 
constraints becomes, in general, impracticable. To overcome 
this difficulty the gradient of these functions will be 
used to predict their values at a new design point and. its 
validity will be investigated. 



(vi) Lastly, to ’evaluate the gradients of the above 
mentioned functions and the gradient of the minimizing 
fuiction at a design point, it tatIII be assumed that for 
small variations in the components of the design vector 
the maximum values of these functions occur in the vicinity 
of the positions at which they have occurred at the base 
point. The validity of this, assumption will be determined. 



CHAPTER II 


KIHSTOBI/ASTODYNAMIC AHAIYSIS 

2.1 Introduction 

When a mechanism runs at a certain speed the resulting 
rigid body inertia forces and the external forces acting on 
the mechanism cause elastic deflections in the links which 
become substantial when the running speed is high. Since the 
rigid body inertia forces as well as the effect of the 
external forces on individual link vary during the motion of 
the mechanism, they produce forced vibrations in the links 
of the mechanism causing the links to deflect further about 
their rigid body position at any instant. The present 
chapter is aimed to make this kinetoelastodynamic analysis 
of the mechanism and determine the stresses in the links 
due to elastic deformations. The elastic accelerations of 
the links obtained during this analysis are then used to 
determine vibratory inertia forces. 

Section 2.2 presents the various assumptions made 
for this analysis and the mathematical model adopted for the 
mechanism. Section 2.3 deals with the element analysis. 

In section 2.4 a procedure for the assembly of element 
matrices into system matrices is presented. The eCLuation 



of motion and a method of determining transformation matrix 
are given in section 2.5 and 2-6 respectively. Section 2.7 
describes a method to eliminate the rigid body degree of 
freedom. Section 2.8 gives the solution of the equation of 
motion and the determination of stresses in the links. In 
section 2.9 the determination of quasi-static stresses is 
presented. Lastly, section 2.10 <^eals with an alternative 
procedure of determining the deflections and stresses by 
suppressing the rigid body degree of freedom. 

2.2 Assumptions and Mathematical Model 

All the links of the mechanism are considered to be 
completely elastic. The important assumptions made during 
the analysis that follow are; 

(i) The elastic deflections are independent of the rigid 
body motion, 

(ii) the deflections are so small that linear theory can 
be used, 

(iii) the load vector, the stiffness matrix and the mass 
matrix remain constant during the small intervals 
into which the complete cycle of motion of the 
mechanism is discretised, 

(iv) the effects of the friction, tolerance, clearances 
and impact are neglected. 

A planar four bar mechanism having straight links with 
uniform cross-sections is considered for presenting the 
procedure of analysis . 



initial rigid' body configura- 
tion of the mechanism at the instant x = The load 

vector, the mass matrix and the stiffness matrix are deter- 
mined at this position which are assumed, to remain constant 
during a chosen interval Ax. Under the action of rigid 
body inertia forces produced by gross rigid body accelerations 
and the external forces the mechanism gets deflected as an 
elastic structural system to a new position O^A 2 B 202 shown 
in Pig. 2.1. Using the displacements and the accelerations 
the stresses in the links and vibratory inertia forces are 
determined, respectively. Prom the deflected position 
^±^2^2^ 2 mechanism is allowed to move as a rigid link 

mechanism through the interval A x to reach the position 
0-,A^B,0-, which is treated as the instantaneous structure 
for the next interval of time Ax. The displacements and 
velocities of the elements obtained at the end of one 
interval are used as the ’initial conditions’ for the next 
interval. Starting from some initial configuration the 
above procedure is repeated until the cycle is complete or 
the steady state condition is reached. 

2.3 Element Analysis 

The links of the mechanism are considered as beam 
elements with six coordinates, three at each end. Out of 
these three, two represent linear displacements, one along 
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the element axis and* the other perpendicular to it. The 
third represents the rotation at the end. The displacements 
at these coordinates are taken as to 5^(1:). To 

represent the continuous deflection of the element (Pig. 

2.2) in terms of these nodal displacements the following 
polynomial relations are assumed 
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where is the continuous axial displacement due to axial 

forces only and u^ is the continuous transverse displacement 

<y 

due to bending only. to are functions of time t only. 
These polynomials agree to the accuracy of the engineering 
theory of static bending of beams and also permit rigid body 
motions in addition to elastic displacements [43]. 

The additional transverse deflection due to shear is 
determined as follows: 


and 
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where, = the transverse deflection due to shear effects only! 

iy ( 

Sj, = the shear force acting at a position x of the link | 
Ag = the effective shear area of the position x 
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G- = the shear modulus of elasticity 

B = the Young's modulus of elasticity 

I = second moment of area of the cross-section. 

Using eqs. (2.3) and (2,4) Sj, is eliminated to get 
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Integrating this equation the substituting equation (2.2) 
into it yields the following 
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The total displacement is given by 
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where the constant of integration C in eq. (2.6) is absorbed 
with constant into in eq. (2.7). 

To determine to Gg the following boundary conditions 
are employed 
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At X = o 

At X = 1 

The resulting transverse deflection can thus be expressed as 
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element . 


Using the engineering bending theory the longitudinal 
displacement due to bending is determined as follows 
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where 5 - y/ 1 ; y = distance of the point from neutral axis. 

and C 2 in eq.uation (2.1) are obtained by applying the 
boundary conditions on axial displacements, i.e. u| = 
at X = 0 and u^ = u^ at x = 1 . The total axial displacement 
u is given by 

(x,y,T;) = uJ(x,t) + u^(x,y,T) ( 2 . 11 ) 

Therefore, the continuous displacement u = { u ,u 1 is 

'• X y/ 

obtained from the nodal displacement u = {u^,...,Ugj 
as follows 
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u = 

au 

(2.12a) 

and as 

such 
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and 
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where the parameter a is called the shape function and is 
given hy 
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According to linear theory the normal strain and the 
sharing strain at any point of the link are obtained in terms 
of deflections as follows 
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Using eq. (2.12) along with eqs. (2.14) and (2.15) the conti- 
nuous element strains at any point of the element in terms 
of the nodal displacement u are given by 
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where the element strain matrix b is 
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Prom the element strains, the element stresses at any 
point of the element are obtained as follows 
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and A = nominal area of cross-section. 

Ihe inertia forces due to rigid body accelerations are 
regarded as the continuous body forces B(? ,x) acting over the 
element. I'he components of these body forces at any point 
of the element are given by 


B(? ,t) = - 
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where p = density of the material, a^^ and a^^ are the 
absolute rigid • body normal accelerations and ^By 



■the tangential accelerations at the end points A and B of 
the element respectively. These accelerations are obtained 
by the kinematic analysis of the mechanisms which is presented 
in Appendix A. 

The element stiffness matrix k, the elemient mass matrix 
m and the element load vector p in the element oriented 
coordinates are obtained using the principle of virtual work 
[42] as follo-ws 


k = 

. , t 

J b X t dv 

(2.21) 


y 


m = 

"t 

/pa a dv 

(2.22) 


V 


P = 

/ a"^ B dv 

(2.23) 


V 


where t indicates the transpose of a matrix. Integrations 
are taken over the whole volume v of the element and dv is 
elemental volume of the element. 

The explicit expressions for these matrices are given 
in Appendix B. 

Since the orientations of various links in a mechanism 
are different it is necessary to define a common coordinate 
system called as global coordinate system (or simply sys-tem 
coordinates). The element matrices are transformed into system 
coordinates from their respective element oriented coordinates 
before they are assembled into system mat ri ces . To achieve 
this transformation, a rotation matrix R is defined as below 
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i 0 
L 


0 




(2. 24) 


where is the rotation matrix for one end of the eleiBent. 
If 0 is the angle made hy the x-axis of the element with 
the X-axis of the system coordinates the rotation matrix R, 
is given by 


R 


1 


cos9 

-sirR 

0 


sine 

cosG 

0 


0 

0 

1 


(2.25) 


I'he nodal displacements u of the element in the system co- 
ordinates are obtained using the following coordinate trans- 
formation 


u = R u ( 2. 26) 

Using the contragradi ent law of transformation, the 
element load vector p, the stiffness matrix k and the mass 

/N 

matrix m in the system coordinates are given by 


p = 

r"^ 

P 


(2-27) 

A 

k = 

R^ 

k 

R 

(2. 28) 


t 




m = 

R 

m 


(2. 29) 


The initial element displacement and velocity in system 
coordinates are obtained in a similar way from the initial 
element displacement and velocity in the element oriented 
coordinates as follows 



3B 


% = % (2.30) 

% . (2.31) 

The above analysis is applied to all the elements present in 
the mechanism. 

2.4 Assembly of Element Matrices 

Many methods are available to assemble the element 
matrices into system matrices for the complete mechanism. 
Some of them are connectivity matrix [28,24j, permutation 
vector method [43,27] and 'code-system' method [44,35]., The 
'code-system' method has the advantage over the others that 
no regular order is to be followed during numbering the 
element coordinates and the system coordinates at any joint. 
Thus any set of coordinates can be grouped together as the 
last set among the generalized system coordinates' which is 
advantageous in placing the coordinates corresponding to the 
rigid body degree of freedom and the fixed ends supports as 
the last set among the generalized coordinates. 

In this method the elements of the mechanism are 
numbered according to any choice. The system, coordinates 
are also assigned numbers as required. An integer matrix 
NS is now defined such that in each row the columns contains 
the number corresponding to system coordinates for that 
element whose identification number matches with the row 
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number in the order in which the element coordinates are 
defined. I'he integer matrix NS for the mechanism shown in 
Pig. 2.3 is given in iable 2.1. 

A sample programme to perform the required assembly of 
the element matrices into the system matrices is shown in 
Table 2.2. 

2.5 Equation of Motion 

Having determined the system mass matrix M, the stiff- 
ness matrix K and the load vector P after the above assembly, 
the special features of the mechanism are incorporated. Por 
example, a flywheel is present at output end of the rocker 
its moment of inertia is added to the diagonal element of 
the mass matrix M corresponding to the system coordinate 
associated with the rotation of that rocker end. Similarly 
the forces or torques exerted at the input end shaft are 
included in the corresponding elements of the load vector P. 

After all the special features of the mechanism are 
accounted for the final equation of motion for the complete 
mechanism is given by . 

MU + ctJ + KU = P ( 2.32) 

for < X < Tq + At, with the initial conditions: 

U(t^) = and UCtq) = Uq (2.33) 
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i'ABLE 2.2. S^lPLE PROGRAM POR ASSEMBLY OP 
ELEMENT MATRICES 

DO 1 I LINK = 1,ME'E 
DO 1 1=1, NEC 
IROW = NS{ILINK,I) 

- SP(IR0¥) = SP(IR0¥) + EP(I,ILI1TK) 

SD(IR0¥) = BD(I,IIINE) 

SY(IR0¥) = BV(I,ILINZ) 

DO 1 J=1,NEC 
ICOL = NS(ILINK,J) 

SM(IR0¥,IC0L) = SM(IR0¥,IC0L)+BM(I, J,ILINK) 
SK(IR0¥,IC0L) = SK(IED¥,IC0L) + BK(I , J, ILIHZ) 

1 CONTI NUB 

where MNB = number of links in the mechanism; NEC = number 
of element coordinates; BP, ED, BV, EM and EK are the element 
load vectors, element initial displacement vectors, element 
initial velocity vectors, element mass matrices and element 
stiffness matrices in system coordinates of the elements of 
the mechanism; SP, SD,SV,SM and SIC are the system load 
vector, system initial displacement vector, system initial 
velocity vector, system mass matrix and system stiffness 
matrix respectively. 
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where 0, U(t) , and 13'^ are the damping matrix, displace- 
ment vector, initial displacement vector and initial velocity 
vector, respectively. 


2.6 Determination of Transformation Matrix 

The equation of motion (2.32) will be solved using 
normal mode method. For this purpose, it is first required 
to determine the solution of the following eigenvalue problem 


Ml + KX = 0 (2.34) 

In the case of a mechanism, the stiffness matrix K 
is singular [28] and the degree of singularity of the matrix 
K is sarnie as the number of rigid body degrees of freedom of 
the mechanism, i'o overcome this difficulty in solving eq,. 
(2.34) the equation is m.odified using the following pro- 
cedure [35]. 


The total set of the coordinates X of eq. (2.34-) is 
divided into a subset called 'slave coordinates' so that 
no rigid body displacement is possible through the remaining 
subset X^ called 'master coordinates'. The total coordinates 
in Xg are Same as the number of rigid body degrees of 
freedom in the mechanism. There are many possible combina- 
tions to form the set of coordinates Zg, but it is more 
advantageous to take the coordinates corresponding to input 
motion as the slave coordinates . The eq. (2.54) is 
partitioned according to these subsets as shown below, 
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M 


•AA 


I 


M 


AB 




A 






K 


AA 


K 


AB 


K. 


BA 


K 


BB 


1 

^ y 'l 

i 

hi 




2.35) 


When a sniall rigid body displacement dXg is given at the 
slave coordinates then the corresponding rigid body dis- 
placements ax the master coordinates are given by the 

relation 


dX^ = 1' dX^- (2.36) 

where i‘ is called the displacement transformation matrix 
[42] and is a function of the geomexry of the mechanism only. 
I'he displacement transformation matrix 1', for any mechanism., 
Can be determined as given below. 


Since no elastic forces are developed when purely 
rigid body displacements take place in a mechanism mathema- 
tically one can say 


K 


AA 


K 


AB 


K. 


BA 


K. 


BB 


dX 


A 


dX- 


(2.37) 


B 


from which one gets 


^AA ‘^•^A ^AB '^^B 


%A '^^A ^ %B d^B 


Substituting e.q.. (2.36) in ,eq 
tions are obtained. 



(2.38) 


= 0 

(2.38) the following rela- 
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= -Km 
Kba'-*-' Kbb 

Eq. (2.39) is used to 
singular and thus its 


^AB 

0 

determine T. ihe matrix 
inverse exists. 


(2.39) 

(2.40) 
is non- 


2.7 Elimination of Rigid Body Degrees of Freedom 


Let a Small virtual rigid body displacements 
given at the slave coordinates. l‘he corresponding 
rigid body displacements at the master coordinates 
by (using eq. (2.36)), 

= T 6Xj^ 


6X3 be 
virt ual 
are given 

(2.41) 


Applying the principle of virtual work to the equili- 
brium eq. (2.35) and remembering 'ohe fact that the elastic 
forces do not produce any work due to rigid body displace- 
ments and work is done by the inertia forces only, the 
following equation is obtained. 


b^A ^ ■*" ^ ^"^BA^A - b (2.42) 




Substituting eq. (2.41) in eq. (2.42) and remembering the 

■ " 2 

fact that 6Xg is arbitrary, and X = -o) X, the following 

relation between X. -and X-r, is obtained 

A h 


X-n = X 


"B 


■A 


where 'i'“ 


. AB ^ ^AA BA^ 


(2.43) 

(2.44) 


^5 


Resolving for small values into two parts [42] as given 
in eq. (2-45) and substituting it in eq. (2.43) the eq. 
(2.46) is obtained 


^ ^Ae (2.45) 

where represents the displacement at the master coordi- 
nates relative to the rigid body configuration when Xj^=0 . 

Xg = (I - I'^T)"^ (2.46) 

I'he size of the matrix (I-T'"’!) is equal to the total number 
of the rigid body degrees of freedom. 


Eliminating X^ and Xg with the help of eqs. (.2.45) and 
(2.46) and making use of eq. (2.39)? eq. (2.35) can be 
written in terms of as follows, 

(-03 2m' + K^^) X^^ = 0 (2.47) 

where w = natural frequency of vibration of the system, 

M' = effective mass matrix in the fr ee-vibration of 
the mechanism and is given by 


M’ = + (M,/i' + M^t^)(I-T-T) V 


A A 


BB' 


2.48) 


The size of the matrices and M' are (n-w)x(n-w) whore 
n is number of system coordinates and w is the number of 
degrees of freedom of the mechanism. In case of single 
degree of freedom system (l-i'*I!) is scalar and therefore Id’ 


is symmetric. Since is non-singular the equation (2.47 
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can. now be solved [45] and the natural frequencies and 
natural modes can be obtained. Let <p ' be the modal matx’ix 
obtained by solving eq. (2.47), the modal matrix cp corres- 
ponding to eq. (2.35) is obtained using eqs . (2.45) and (2.46) 
as ' given below 



/ ^ ] 

! 9 b J 


where cp-g = (l-I''"T) ^ (2.50; 

and 9 j\ = 9 ' + 4*2 (2.51; 

■2.8 Solution of Equations and Determination of Stresses 

I'he vector U in eq. (2.32) allows the rigid body 
motions along with the elastic motions (by virtue of the 
polynomials in eqs. (2.1) and (2.2 )), but the modal matrix 9 
contains only the elastic modes and no rigid body modes. 

Thus, to determine the solution of the . equation of motion , 
(2.32) by normal mode method the system coordinates U are 
transformed to normal coordinates with the help of the 
following relation 

U = cpT) : ■“ :( 2.52a) 

where 11(1:) represent the displacement vector in the normal 
coordinate system. ; 

Since 9 is assumed to remain constant during the chosenj 
interval At, the following, relations may readily be obtained 

^Abstract of the method from this reference is presented 
’ in Appendix 0. 


M 


for velocity and acceleration in normal coordinate system 
by successive differentiation of eq. ( 2 . 52 a) 

(2.52b) 

and U = 91) (2.5 2e) 


Plugging eqs . ( 2 . 52 ) into eq. ( 2 . 3 , 2 ) and, premultiplying 
ty 9^ y the following equation results 


where s 


Mp + Of] + En = P 
M = cp Pkp 5 C = 9 C9 

K = 9^ K9} p = 9^P 


(2.53) 


(2.54) 


Using eqs. ( 2 . 53 ), the initial conditions are given by 


% = (9^9)"^ 9"^ . (2-55) 

% = (9^9)”^ 9^ (2.56) 


Since the modes in 9 are orthogonal with respect to 
the matrices M and K, M and K are diagonal matrices. If 
daping matrix C is taken proportional to either Pi or K or a 
linear combination of them, the matrix C also becomes a 
diagonal matrix, and eq. (2.53) is decoupled. I'he standard 
practice is to take C as proportional to the critical 
damping of the system at the normal modes [ 42 ] . Ihus the 
following uncoupled equations are obtained 


m. n. + 2£. m. 
1 i ‘"i 


Hi .+ k. q. 


Ps 


(2.57) 


i — 1 , 2 j 


, m. 
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where , 

m 

= 

number of normal coordinates 


^i 

= 

displacement at the i-th normal coordinate 



= 

velocity at the i-th normal coordinate 



=: 

acceleration at the i-th normal coordinate 


m. 

i 


i-th diagonal element of Fi 


k. 

X 

= 

i-th diagonal element of ii 


Pi 


i-th element of P 




i-th natural frequency 

2m^ u 

i 

= 

critical damping at the i-th normal mode 



= 

damping ratio (actual damping to critical 

damping) at the i-th normal mode. 


In general the nature of system taping is such that 
the effect of higher modes is suppressed, therefore, the 
common practice is to take higher values of for the 
successively higher modes. The damping matrix C as taken 
here in the eq. (2.53) can be expressed as follows 


c = 2 


“ 1^1 


[M] (2.58) 


j m [ 

Since it is assumed that load vector P (and as such P) . 
remains constant during the interval , the solution of 

eq. (2.53) is obtained as shown below 

-'3 ■ s ^i / 

p.(r) = e ^ (I. cos p-r' + J.sin p. t ' ) + g- (2.59) 

1 1 X X i 



-19 


ri^(T) = e 


and 


-a . T ' 
1 


( 2,b0) 


P- p 

“ ~ Ti- - 2 « .C . -n. 
Dlj_ 1 1 1 


(2.61) 


i = 1, 2, . . , 


for X < T < T + . 


where x 


X - x^ = ^x 


O’ . =: f . 03 • 

1 11 


lii = n 


h 


n° - PiAi 

1 = ('Oi + 


n° = i-th element of the initial vector t]^ in eq. (2.55) 
and q? = i-th element of the initial vector in eq. (2.5b) 

The displacement q, velocity q and. acceleration x] .as 
obtained from eqs . (2.59) to (2.61) are transformed into sys- 
tem coordinates using eq. (2.52). i'he element displacement, 
velocity and acceleration vectors in system coordinates are 
obtained from, the system vectors using integer matrix hS. 
'these vectors are then finally transformed into element 
oriented coordinates with the help of eq. (2.26) . The dis- 
placements thus obtained are used, to determine the stresses 
in the links of the mechanism with the help of eq. (2.18). 

The acceleration vector obtained is Used to determine the 
vibratory inertia forces. 
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2.9 Determination of Quasi-Static Stresses 


If the load vector corresponding to the master coor- 
dinates calculated at any instant is allowed to act on 
the mechanism as static forces the elastic deflection 
produced in the mechanism in the system coordinates (assuming 
the mechanism fixed at the input end. i.e. the crank as 
cantilever) is given by 


U 


static 


= K 


-1 


AA 


■A 


( 2 . 62 ) 


This displacement vector is transformed into element 
displacement vector in -system coordinates using integer 
matrix IfS . The element displacement vector thus obtained 
is transformed into element oriented displacement vector 
using eq.. (2.2b). Plugging these element displacements in 
eq. (2.18) .the quasi-static stresses in the links of the 
mechanism are obtained. 


2.10 System Analysis by Suppressing Eigid Body Degree of 

Freedom 

It has been pointed out in section 2.6 that the matrix 
k in eq. (2.34) is singular. Thus to solve eq. (2.34), a 
way of modifying che equation by eliminating the rigid body 
degree of freedom has been presented in the previous sections 
Another alternative is to treat the input crank end. as 
fixed, like a cantilever beam, thereby suppressing the rigid 
body degree of freedom. Tills is achieved by deleting the 
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rows and columns of the matrices M and K corresponding to 
the coordinate representing rigid body degree of freedom 
and' also by deleting the elero.ent corresponding to this from 
the vector U in equation of motion (2.32). I'hus, the 
resulting stiffness matrix denoted by is non-singular 

.Hjn. 

and the equation of motion can now be written as 


^ ^AA^A ^AA ^^A ~ ^A 


(2.65) 


where is the resulting mass matrix, the damping 

matrix for this problem, the corresponding load vector 
and U the resulting displacement vector. 

A 


Bq. (2.6 3) can novf be solved using normal mode m,ethod. 
Bor this purpose solution of the following eigenvalue 
problem is sought 


^^AA ^A ^AA ^A 


0 


(2.64) 


Since is non-singular the above equation is 

solved as before and the natural frequencies and natirral 


modes are obtained. The rest of the procedure is same as 
outlined in section 2.8, and the stresses in the links of 
the mechanisms are obtained. 


hi 


CHAP'JJER III 


KINEMA'IIC SWHESIS AED DYNAMIC EORCE ANALYSIS 


3.1 Introduction 

A method of kinematic synthesis for a planar four bar 
mechanism is presented in section 3.2. 5Co determine the 
maximum unbalance shaking force and shaking moment during 
a complete cycle of motion of the mechanism and also to 
determine the amount of the input torque required to maintain 
equilibrium at each interval into which the complete cycle 
of motion of the mechanism is discretised, it is required 
to make complete dynamic force analysis of the mechanism at 
the configuration corresponding to each intervai. iSection 
3*3 deals with a procedure for this dynamic force analysis. 

In section 3.4 a method is developed to make the dynamic 
force analysis when counterweights are attached to the 
links in order to redistribute their masses with the aim 
of reducing the shaking force and/or the shaking moment. 

In section 3. 5» a procedure is developed to determine the 
contribution due to the KED inertia towards the shaking 
force and shaking moment. Section 3.6 deals with the 
determination of system mass matrix when counterweights 
are attached at the ends of the links. 
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3.2 Kinematic bynthesis 

I'he kinematic synthesis of mechanism has received a 
great deal of attention ill to date. A number of methods 
for the kinematic sjmthesis are now available [46-48] . The 
present study does not aim at discussing the merits and 
demerits of these methods, bince the complex number method 
of the kinematic synthesis is most widely used for planar 
mechanisms, it is adopted in the present study. 

Though the method is a general one and can be used 
to synthesise a planar mechanism for any use, i.c. path 
generation, function generation or rigid body guidance, it 
is discussed here with reference to rigid body guidance. 

buppose it is proposed to guide a rigid body through 
finitely separated positions. Let the firm lines in Pig; 3. 
show the intial configuration when the rigid body is in its 
first position at coincident with the coupler point P. 

and 0^ are the fixed supports and 0 is any base point. 

Zq to represent the vectorial positions of the first 
support from the base point and the various links of the 
mechanism as indicated in Pig. 3.1. The vector Zg repre- 
sents the position of the coupler point p from the second 
joint of the coupler, shown by dotted line in Pig. 3.1. 
let the vectors P.^ and R^ represent the three finitely 

separated positions of the rigid body attached to the 
coupler point P. 




The loop closure equations for the point P at the 


initial position is written in the following two ways 


Zq + 22 + 


= R- 


1 


(3.1) 


^0 + ^1 ^4 ^6 = h 


When the mechanisin moves and the rigid body occupies 


the jth position, let the links be rotated through 0 . ( cranM 


. (coupler) and 9. (follower) from their initial positions 
J J 


respectively, then the above loop closure equations for the 

if. 


jth position are 
'o ■ -'2' 


Z„ + Z^e ^ + Zc-e ^ = R. 




2 o + ^1 + Z^e 


i 

^ + Z^e ^ 


( 3 . 2 ) 


R. 


Subtracting equations (3.1) from equation (3.2) the follOwi]| 
relations are obtained. 


ie . iif.. 

Z^ (e J - 1) + Z. (e - 1) = r, 


i9, . iif.. 


(3.3) 


Z. (e ^ - 1) + Z. (e J - 1) = r . 


where r. is the relative displacement vector of the j 
J 


position of the rigid body from its initial position . For 
three rigid body positions j is equal to 2 thus one has 


i0-) 

Z2 (e ^ - 1) + (e ^ ~ 1) 

Z2 (e ^ 1) + 2^ (e 2-1) 

i tPi i 'I't 

ie - 1) + (e ^ - 1) 

iTj)p 

Z4 (e 2 _ 1) ^ 2g (e 2 _ 1) 

I'hase are four equations in four complex variables. I'he 
ao-gular rotations of the links may he arbitrarily chosen or 
sone of them may be specified in the problem, generaily 0's 
are specified, and Z25 Z^, Z^ and Z^ are deteirniiHfed Using 
oqviaticns ( 3 . 4 ). and are then determined using the 

links closure equation and the closure equation of the 
couple:- triangle as given below 

Z? Z3 = Z-i^ 

( 3 . 5 ) 

+ 2^ = Z^ 

Therefore the equations ( 3 . 4 ) and ( 3 . 5 ) yield all 
the link lengths of the mechanism to serve the required 
aim. However, the number of possible solution are infinite 
since arbitrary values can be assigned to i|; * s and tp's . 

3.3 Dynamic Force Analysis 

At each stage of iteration during the optimization 
process it is required to determine the maximum shaking 
force, the maximum shaking moment etc. during a complete 


= 


= r 


1 


5G 


( 3 . 4 ) 
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cycle of motion of the mechanism. Ihe following procedure 
is adopted for this purpose. 

In Chapter II a transformation matrix has been 
obtained, (section 2.6) which relates the small rigid body 
displacements at the 'slave' set of coordinates to the 
corresponding small displacement at the 'master' set of 
coordinates. The sam.e logic is extended when the coordinates 
representing the two translational rigid body degrees of 
f^i-^edom, at each fixed, hearing are also included [35 J- These 
four cO “Ordinat es , two at each fixed bearing, that allow 
additicnal rigid body degrees of freedom, are grouped 
togsuher with the rotational degree of freedom at the input 
cranlc end, placing them as the last system coordinates. 

This set of coordi.nates forms the group of ''slave' cocr'dl- 
nates and is represented by the subscript C. The remaining 
set of coordina.t ss , xrliicli are termed as 'master' coordinates, 
is denoted by the subscript A- 

At any given configuration the mechanism is in equili- 
brium under the action of the external forces acting on it, 
the inertia forces, the torque applied at the input crank 
end and the reactions at the fixed bearings. Applying the 
principle of virtual work and remembering the fact that the 
elastic forces do not produce any work during, rigid body 
displacements the following equation is readily 0bta.ined 

^0 “ ( 3 - 6 ) 
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where Tp is the vector of the unknown, reactions at 
coordinates, is the known force vector that include the 
inertia forces and external forces at the slave coordinates 
and is the known inertia forces and the external forces 
at the master coordinates. The transformation matri ^ 
which relates the small rigid body displacements at the 
slave coordinates to the corresponding small 
displa.cement s at the master coordinates i® st a' 

• , .T-j-in. “sections. 6, 

ing a procedure similar to that described, m 


giving as follows 


T 


R 


"^AA ^AO 


(3-7) 


where and are the submatrices of the par 

stiffness matri 


: ri X TT 


[K] 


K. 


'AA 

AC 

CA 

^CC 


(3.8) 


. eQ • ' ( 3 • ^ ) t he 

Having determined the vector Tq ^ 

+• a noint midway 

shaking force and the shaking moment aboui a r 
between the two fixed supports are .given by 

(3.9) 


PS 

= 

MS 

1 

= _£ ( i£i 

2 ^ ( 



'c5 ‘^c3'^ 

where PS is the shaking force, MS 


(3.10) 

the net shaking moment 


exerted on the foundations, 1 the distance 

’ o 


Petween the two 
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attached to the link respectively. Let the suffix A and B 
stand for the parameters of the counterweights at the left 
and right ends of the link respectively. Let o) and a are 
the angular velocity and angular acceleration and a. 

and a-Q^j are the linear accelerations in the x and y 

directions at the left and right ends of the link respec- 
tively. 'i'he additional nodal inertia forces due to the 
counterweights are then given by 


Pi = 


^"'AX 

r 

A 


P2 = 


^^Ay - 

^A 

a) 

P5 = 


- P2 

r 

A 


P4 = 


^^Bx ■ 

^B 

-2) 

P5 = 

-®B 

(a^y + 

^B 

a) 

Pb = 

-IbQ 

: + P5 

^B 



where P2_ P6 component of the additional nodal 

inertia forces due to the counterweights acting at the 
six coordinates of the element (link). Ihese nodal forces 
are added to the previously ' calculated element nodal forces 
due to the inertia of the link alone. This resultant set of 
the nodal forces gives the element load vector in element 
oriented coordinates due to the combined rigid body inertia 
of the link and its counterweights. 
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fixed supports, is the first element of the vector, I'p, 

which corresponds to the input torque at the crank end and 
^'c2 ^'c5 remaining elements of the vector 1‘ 

u 

which correspond to bearing reactions. 

3.4 Dynamic Force Analysis with Counterweights 

When the masses of the links are redistributed by 
placing counterweights in order to reduce the shaking force 
or shaking moment, the inertia forces acting at the nodes, 
determined using the finite element approach, due to the 
links of uniform 'cross-section, are modified. Ihe addi- 
tional nodal forces due to the inertia of the attached 
counterweights are determined as follows. 

I 

The counterweights that are attached to a link are 
generally either circular with circumference passing 
■ through the center of the joint of the link, or the sector 
of a circle with its centre coincident with the center of 
the joint or simply an extension of the link beyond the 
joint having similar cross-section as the link. Figure 3.2 
shows a link with counterweight attached at both ends. A 
link may not have counterweights at its both ends but to 
generalize the presentation counterweights are assumed to 
be present at both the. ends. Let m, r a.nd I represent the 
mass, the distance of the centroid from the joint of the 
link and the mass moment of inertia of the counterweight 



G1 


section 3.4 for determining the contribution towards the 
rigid body inertia forces due to the counterweights. The 
only change made in the equations (3.11) is that the 
accelerations due to elastic motion are substituted in 
place of the rigid body accelerations. Thus one gets 


I — 1 
© 

p. ■ 

= 

“a 


^A 

^e3 

Pe2 

= 

““a 


^A 

%3 

Pe3 

= 

-^A 


- Pe2 

^A 

Pe4 

= 

-mB 

‘**64 

^B 

;^e6 

Pe5 
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where j 

of the KED inertia load vector due to counterweights, 
velocity vector v and acceleration vector a of the element 
(link) in the element oriented coordinates, respectively. 

The other variables are same as explained in section 3*4. 

The velocity vector and acceleration vector mentioned above 
are the ones that are obtained at the end of the previous 
interval by the kinet oelastodynamic analysis of the 
mechanism. 


3.5.3 Determination of Total lED Inertia Do rces 

Having detenained the KBD inertia load, vector as 
above it is added to the element resultant load vector 
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determined in section 3.4. This final element inertia load 
vector in element oriented coordinate is transformed into 
element load vector in syst.em coordinates using the eq,. (2.26). 
These transformed elements load vectors are assembled to 
yield system inertia load vector. The inertia load vector 
thus obtained is added to the system inertia load vector 
determined in section 3-5.1 to give the final system inertia 
load vector due to the combined rigid body and KBD inertia 
for the complete me chanism, • i . e . including all the links 
and their counterweights. It is this inertia load vector 
which is used to determine the shaking force, the shaking 
moment and the input end torque at any configuration of the 
mechanism as outlined in section 3.3* 

3.6 Determination of System Mass Matrix with Counterweights 

When the counterweights are attached at the ends of 
the links the system mass matrix M in eq. (2-33), obtained 
considering the regular links with uniform cross-section, 
needs modification. To achieve this it is assumed that as 

if the whole of the counterweight mass is situated at- the 

; 

joint of the link near to that counterweight to account 
its contribution corresponding to the translational coordi- 
nates of the joint. To account the contribution of the 
counterweight corresponding to the rotational coordinate at 
the joint the moment of inertia of the counterweight about 
the center of the joint is determined. These masses or 



moment of inertia are ad^ed to the corresponding diagonal 
elements of the element mass matrix in element coordinates. 
The element mass matrices, thus obtained, are first trans- 
formed into element mass matrices in system coordinates, 
and they are then assembled to give system mass matrix M. 



CHAPTER IV 


optimzation problem ahl solution procedure 


4.1 Introduction 

In mechanism design one of the two methods may be 
followed, namely; (a) Direct design methods, (b) Mathematical 
programming techniques. In the former procedure a set of 
equations is formed based on design criteria (kinematic or 
dynamic or both) relating the design variables. These 
equation are then solved by various means to find the design 
variables. This procedure is not suitable if the required 
design has to satisfy certain inequality constraints. One 
alternative is to design a mechanism using this procedure 
and see whether the constraints are satisfied. If they 
are violated the design is revised by trial and error. 
Eurthermore for complex pro bl en© , where there are a number 
of possible design objectives, including both kinematic 
and dynamic performance criteria and where there are a 
number of design constraints this procedure is not applic- 
able . In such situations it is advantageous to adopt 
mathematical programming techniques. Using this technique 
it is also possible to make a trade-off study between 



various competitive objectives. The application of mathe- 
matical programming techniques to problems in mechanism 
design has been presented in a review paper by Fox and 
Grupta [49] • I'hey have discussed the available tools of 
optimization as applied to mechanism design. In the 
present study, where two stages of mechanism design, namely, 
balancing of mechanism and strength considerations are 
both combined into a single design objective, this is the 
most convenient procedure to be followed. Therefore the 
problems in the present study are solved as nonlinear 
programming problems. 

4.2 nonlinear Programming Problem 

A nonlinear programming problem can mat hematically be 
posed as; 

Mnimize, 

f(X) 

s ub j e ct 1 0 , 

S (x )<0 i = l, 2 ,...,m 

and , 

Ij.(X)=: 0 j = l, 2, ..., p 

where f is the objective function of n variables X to be 
minimized, g. are m inequality constraints and 1. are p 
equality constraints which must be satisfied at the optimum 
point . 



4.5 Problems Undertaken Luring Present Study 

As pointed out in Chapter I the present study is 
categorized into two types of problems and hence the follow- 
ing two nonlinear programming problems are formulated . 

4.3.1 first Problem 

In this problem it is required to synthesize a suii- 
able mechanism which serves a specified aim (such as rigid 
body guidance, path generation or function generation) wifh 
optimal dynamic characteristics with reference to unbalance 
in the mechanism. Ihus the objective is to minimize the 
shaking force and/or the shaking moment transmitted to the 
machine foundation. 

The design parameters for this problem are the unknown 
relative angular positions of the coupler and the follower 
when the crank moves through given positions (angles) to 
guide the rigid body through required positions. The 
cross-sectional areas of the links are taken as variables . 

A single proportionality factor is taken for this purpose 
for all the links in order to limit the number of variables. 

The major constraints Imposed on the optimization 
process are as follows; 

(i) from kinematic point of view the mobility 
conditions must be satisfied and also the transmission 
angle at any position of the mechanism should not be less 
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than a specified minimum. Assuming the mechanism to be 
crank rocker type the mobility conditions (knovm as Grashof 
criteria) and other constraints are expressed in terms of 
link lengths and transmission angle X as follows 


(dl+12) “ ^ 

1^-14 - ~l2 < ^ 


I2-II -< 

1^-1 . < 
2 min 



90-X 


max 


- 0 


min 



(4.1) 


where 1^^ (i = 1, 2,5,4) are the lengths of ground link, 
crank, coupler and follower respectively, is the 

minimum permissible length of the crank and 6 is the 
least permissible value of the maximum deviation of the 
transmission angle from 90°. 'i'he first three of the above 
constraints are Grashof criteria, the fourth one takes into 
account the minimum practical limit for the crank length 
and the last - one is from the transmissibility consideration 
of the forces. 


(ii) The stresses in the links are kept within a 
specified maximum limit to safeguard against possible 
failure due to want of strength. Mathematically it can 
be expressed as 



a < 

max 


0 


(4.2), 


i = 1,2,3 



where are the maximum stresses during a complete cycle 
of motion of the mechanism produced in the crank, coupler 
and follower, respectively, is the maximum .allowable 

stress for the material of the links. 

,(iii) The deflections of the strategic, points in 
the mechanism are not allowed to exceed beyond a specified 
limit so that the required goal is met- within 'a tolerable 
accuracy. However, in the solution of present problem no 
constraints are imposed on deflections. 

An example has been solved taking into account the 
effect of KIED inertia during the determination of the shaking 
force and the shaking moment transmitted to the machine 
foundation. In another example the contribution due to 
KED inertia is not included and it is solved taking into 
consideration only the contribution due to rigid body 
inertia to investigate the effect of KBD inertia on the 
design. 

4.3.2 Second Problem 

For a given mechanism '( i . e . with fixed link lengths) 
the link masses are redistributed by placing counterweights 
to achieve optimum dynamic characteristics with reference 
to unbalance in the mechanism- Therefore, the objective in 
this case also is to minimize the shaking force and/or 
shaking moment transmitted to -the foundation. 



The design variables are the size and position of 
counterweights to be attached to the links. The cross- 
sectional areas of the links are allowed to vary. To 
reduce the number of variables, the cross-sectional areas 
of all the links are varied in the same proportion by taking 
a common proportionality factor. 


The major constraints considered are: 


(i) The stresses in the links are constrained to 
take a value less than a specified maximum permissible value 
so that the links have sufficient strength. Mathematically 
it can be represented as 


0 . 

1 


max 


0 


(4.3) 


i = 1, 2, 5 

where a. and a are same as defined in 4.3.1. 

(ii) If it is desired to limit the deflections of 
vital points in the mechanism then constraints can be 
imposed on them. However, in the present problem these 
constraints are not included. 

(iii) In case the shaking force alone is minimized 

a constraint on shaking moment and input torque are imposed. 

In this case also in one example the contribution of 
KBD inertia towards the shaking force and shaking moment 
has been considered. However to estimate the effect of 
KED inertia the same problem is solved excluding KBD inertia 
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Berkof and Lowen [6] showed that in a four har meclia- 
nism with a coupler having symmetry about the line of joints 
the shaking force could be minimized by placing appropriate 
counterweights at 'the end of the crank and follower jointed 
to the fixed supports in such a way that their centroid lie 
on the line of joints of the respective link produced, ibince 
in the present problem a symmetric coupler is taken, the 
centroids of the counterweights attached to the crank and 
follower are kept at 180° with line of joints of the respec- 
tive links. These counterweights are taken as circular with 
their circumference passing through the center of the pin 
of the joint of the respective link. Two sector type weights 
are attached at the coupler ends. The radii as well as the 
angles of the sectors are taken as variables. These counter- 
weights are taken result of an observation in which it 
is noted that they have a favourable effect on coupler 
stresses. As such, the design variables for this problem 
are the radii of the four counterweights and the angles 
of the two sector type counterweights. 

4.4 Procedure Adopted for Evaluating the Haximum Shaking 
Force etc. and Maximum btresses 

The shaking force etc. as well as the stresses in 
the links of a mechanism are not only functions of the 
various parameters of the mechanism but, they also depend 
upon the configuration of the mechanism and as such they 
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continuously vary during the complete cycle of motion of 
the mechanism. It is not possible to express them explicitly 
as, a function of input angle. The general practice is to 
calculate them at discrete positions during the complete 
cycle of motion. The accuracy in the values of the shaking 
force etc. and stresses depends upon the size of the 
interval. A smaller interval size gives a better accuracy 
but requires a larger number of calculations during a cycle 
thereby incurring more computer time in determining the 
maximum shaking force etc. and the maximum stresses in the 
links. Thus a compromise is made between the accuracy 
desired and the available time, furthermore when the 
parameters of the mechanism are changed, not only the 
magnitudes of the maximum shaking force etc. and maximum 
stresses in the links are changed but the positions (i.e. 
angles) of the crank at which these maxima occur also change. 

4.4.1 Determination of Maximum Stresses 

The determination of the maximum stresses in the 
links of a completely elastic mechanism, is a time consuming 
process. Moreover most of the methods available todate 
yield unsteady state stresses initially. To obtain the 
steady state stresses the process has to be repeated for a 
number of cycles which differ in different cases. The 
determination of stresses even for one cycle is very time 
consuming therefore the calculation of steady state stresses 
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by continuing the analysis for many cycles cannot be justi- 
fied in an optimization problem where the minimizing func- 
tion or its gradient has to be evaluated a large number of 
times. Therefore the determination of stresses is restricted 
to first cycle of motion only. Invariably in most of the 
cases the steady state stresses are smaller than the first 
cycle stresses and the design will be safe. 

I'o cut down the computational time further the 
gradients of the stresses in the links are determined at 
the initial point and then while taking various steps in 
the design space during one dimensional minimization to 
find out the optimum step length these gradients are used 
to determine the subsequent values of the mazimum stresses 
in the links. These type of approximations have been 
successfully used in a number of other problems. It is 
found to work well in this case also. 

4.4.2 Determination of klaximum Shaking Force etc. 

For problems in which contribution of KED inertia 
towards shaking force and shaking moment etc. is included, 
the shaking force, the shaking moment, and the input torque 
are determined at each position of the mechanism into which 
the cycle of motion is discretised while determining stresses , 
using the method outlined in Chapter III. The maximum 
values of these functions during the complete cycle of 
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motion are noted. As explained in section 4.4.1 the kineto- 
elastodynamic analysis of the mechanism is a time consuming 
process the gradients of the maximum shaking force, the 
maximum shaking moment and the maximum input torque are 
used, to predict these functions at subsequent design points 
during one dimensional minimization. These approximations 
are also found to work successfully. 

In case the contribution of kED inertia is not taken 
into consideration the shaking force, the shaking moment 
and. the input torque are obtained using direct equations as 
outlined in Appendix E at intervals into which the complete 
motion of the mechanism is divided.. A larger interval may 
be taken for this purpose. Having determined the maximum 
value of the shaking force and the position, i.e. the crank 
angle at which it occurs, a finer interval is taken and 
the shaking force is calculated about this crank position' 
and a more accurate value of the maximum shaking force is 
obtained. If desired this procedure can be adopted to 
obtain accurate values for the maximum shaking moment and/or 
maximum input torque. 

4.5 Solution Procedure 

The procedures available to solve a constrained 
minimization problem fall into two categories: (i) convert- 
ing the problem into an unconstrained problem by suitable 
means, (ii) direct methods [50,51]* 



71 


Without going into a discussion of their merits and 
demerits, which is not the purpose here, the interior 
penalty function method of unconstrained minimization is 
used that fall unider the first category, i'he interior 
penalty function method has an advantage over the exterior 
penalty function method that one, can stop at any level of 
optimization process with the design s.till in the feasible 
region . 

The interior penalty function for problems having 

I 

inequality constraints only is defined as 

cp (X,r) = f(X) - r f — - — (4.4) 

i=l g^(x) 

where r is penalty function ' parameter . 

The function cp is minimized as unconstrained minimiz;a- 
tion problem over X for a fixed positive value of r. In the 
next step r is decreased and 9 is minimized. The’ process 
is continued sequentially till there is no improvement 
during the minimization process. It is reasonable to take 
such values of r initially which makes the function 9 about 
1.5 to 2 times the objective function f. Any suitable 
reduction factor is taken to reduce the value of r for 
the next step. Most commonly a reduction factor of 10 
is taken. 

As pointed out earlier, gradients of the maximum 
stresses in the links shall be used to predict them at the 
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neighbouring design point so as to reduce the computational 
time. Iherefore, it is appropriate to use a gradient method 
for unconstrained minimization. I'he Davidson-Plet cher- 
Powell (DPP) method (also known as variable -metric method) 
has been adopted to solve the problems under study. I'his 
is very effective even for highly eccentric and distorted, 
functions. Golden section method of one dimensional minimiza- 
tion is used, to determine the optimum step length in a certain 
m-Ove direction. I'his method is preferred over the cubic 
interpolation since it req.uires lesser number of overall 
functions evaluations thereby incurring lesser computer time 
in determining optimum step length. 

4.6 Determination of Gradients 

I'he gradients of the maximum shaking forces etc. and 
the maximum stresses are used to determine these functions 
at Bubseq,uent design points. Moreover, the DPP method of 
Optimization requires the determination of the gradient of 
the penalty function at each iteration of optimization. 

Since these functions are not available in the form of 
explicit equations, it is not possible to have closed form 
expressions for these gradients. Therefore, the finite 
difference method, of gradient calculation is resorted to. 

The backward difference formula is used for this purpose. 

If the number of design variables is large sufficient 


computational time is needed when the stresses in the links 
are calculated, each time a design variable is perturbed, upto 
the complete cycle of motion of the mechanism, lo overcome 
this problem the kineto el asto dynamic analysis is ma.de for 
the complete cycle of motion of the mechanism at the base 
point only and subsequently, when each design variable is 
changed by a sm.all amount in turn, the kinetoelastodynamic 
analysis is, made up to a few degrees beyond that crank angle 
at which the maximum stress in the links or the maximum 


shaking force or the maximum shaking moment or the maximum 
input torque (whichever corresponds to largest crank angle) 
has occurred at the ba.se point. This procedure of determin- 
ing the graduents has worked well and has saved sufficient 
computational time. 



CHiPi'BR V 


fiESUITS AMD DISCUSSION 

5.1 Introduction 

The optimization problems formulated in Chapter IV 
are solved by taking numerical examples . Several computer 
programs have been written for tints purpose, which are run 
on D.SC 1090 computer to solve these examples. However, 
before taking up xhe actual optimization problems a study 
has been made to observe the effect of various mechanism 
parameters on the results of the kinetoelastodynamic analysis 
of the mechanism. This analysis is also expected to be use- 
ful in arriving at reasonable values of initial parameters 
of the mechanism for the optimization problems. Therefore, 
the analysis results are presented first in section 5.2 
and then the results for the optimizations are given in 
section 5.5. 

5.2 Analysis Results 

The following four bar 'mechanisms are solved for 
the purpose of studying the effect of various parameters 
of the mechanism.. Any system of units can be chosen for 
the numerical examples. However, DPS system of units 
has been adopted in the present study so that a direct 
comparison can be made 'with the previous works, wherever 
possible . 
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Mechanism A 

Length of the crank = 12.0"; length of the coupler 
= 36.0", length of the follower = 30 .0 • , length of the 
grouna link = 36. u". Ihe other parameters, unless Specified 
otherwise are as follows; material aluminium of sp.^t 
0.101 Ibs/in^ and Young's Modulus = 0.1 x 10® Ihs/in^- 
croBs-section of all the links = 1" x 1"; material damping 
ratio = 0.05 5 angular interval in which the complete 
cycle of motion is discretized = 4°; speed of the input 
crank (assumed constant) = 300 rpm. 


Mechanism B 

Length of the crank = 14 . 42 "; length of the coupler 
= 29.28"; length of the follower = 21.26"; length of the 
ground link = 26.26 ". Other parameters are the same as for 
mechanism a. 


5.2.1 Kinetoelastodynamic Analysis 

Table 5-1 shows the effect of the size of interval 
in which the complete cycle of motion of the mechanism is 
discretized for computational work. The balancing torque 
present at the input end is determined considering the 
rigid body ( RB) inertia forces only. The analysis is 
continued for the first three cycles. A change in the 
^angula,!” interval (A0 ) from 5 to 4 causes an increase in 
the link stresses specially that of the crank during f he 



able 5.1 Effect of Angutar Interval { A 0) 

(% - 0 05 . Links X -sect ion=1X 1 , RPM = 300 ) 
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firsb cycle of motion. 'I'his couia be due to the reason that 
some useful informaiiion might have been suppressed by choos- 
ing a larger interval ( A9 = 5°). However, the difference in 
the link stresses is not much when the interval is further 
reduced, from 4° to 2°. Pigs. 5.1 and 5.2, which give an 
account of the coupler stress versus crank angle, for the 
intervals of 4° and 2° respectively, further enhance this 
conclusion. In examination of the contribution of KED 
inertia towards the shaking force shows a continuous 
increase with the reduction in the interval. It can also 
be noted that the contribution of EDD inertia towards the 
shaking force is not only comparable to the contribution 
of rigid body inertia towards the shaking force but the 
former is even higher than the latter for finer intervals. 
Furthermore 5 the steady state is reached, as can be observed, 
during the second cycle of motion, and the stresses are 
more or less settled to the sam.e maximum magnitude in all 
the three cases. 

I'able 5.2 shows the effect of the inerval sise when 
the contribution of KBD inertia towards the input end torque 
is also included in the analysis. Looking at fables 5.-i- 
and 5-2 it can be concluded that the addition of torque 
at the input end due to KBD inertia, has no significant 
effect on the results. Ihis is because that the contribu- 
tion of KED inertia towards the input end torque is very 
small in comparison to that of rigid body inertia. 




■PointSlressrvsCrank Rotation 
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lable 5.3 gives a comparison between the BSD analysis 

results when diff‘pyc>n+ ^ t ,, 

values for the material damping ratio 

are assumed during the solution process. It is interesting 

no be lohat although there is an increase in the link 

during the first cycle of rotation with a decrease 

in damping ratio, but the maximum crank stress is settled, 

less, to the same steady state value in all the 

three cases considered. However, the coupler and follower 

stresses are settled to slightly higher values with lower 

damping latio. ihe effect of damping ratio on the stress 

distribution in the links can be fully understood by observ-- 

ing ligs . 5,1 which depict the coupler stress 

against the crank angle for damping ratios of 0.05 and 0.01? 

respectively, it can be- observed that the stress fluctuates 

about some mean value. I’he fluctuations are more pronounced 

when bhe damping ratio is smaller. Observing the stress 

during the second, cycle of rotation, it can be concluded, that 

bhe fluctuations are reduced in the steady state as compared 

to the fluctuations during the first cycle of motion which 

depicts the transient state. This beljaviour is common to 

all dynamic response problems. The contribution of IBI) 

inertia towards the maximum shaking force is more for 

lower damping ratio. Furthermore, the maximum shaking force 

in the steady state is more than that in the transient state, 

i<e. first cycle of motion. 
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Table 5-4 shows the results with different assumed 
material damping ratios when the contribution of EID inertia 
towards input end torque is also included. I'he results 
have the oame trend as for the case when the torque due to 
rigid body inertia only is considered to be acting at the 
input end . 

Table 5-5 presents the effect of the crank speed of 
nhe mechanism on the results of kinetoelastodynamic analysis. 
As the running speed increases the stresses in all the links 
are increased more or less in a ratio wliich is the square 
of the speed ratio. Tliis is in accordance with normal 
expectations. Since the acceleration is proportional to 
the squane of the input speed the load vector, which 
consists of forces due to inertia only, is increased as the 
square of the speed ratio. Thus, the increase in the link 
stresses is also of the same order. The shaking force and 
the input end torque due to rigid body inertia are also 
correspondingly increased as the square of the speed ratio. 
The shaking force and the input end torque due to IfED 
inertia are increased in a proportion which is even higher 
than the square of the speed ratio. Here one can observe 
that it is extremely important to take KED inertia into 
account at high speeds. Since the input torque due to 
KED inertia become very high with the increase in speed, the 
link stresses, when the torque at input end due to KED 
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inertia is also included, are much higher than those when 
the torque present at the input end is due only to the 
rigid body (RB) inertia. 

Ihe effect of varying the cross-section of the links 
on the maximum link, stresses, the maximum shaking force, and 
the maximum input torque is presented in Table 5.6. The load 
vector due to rigid body inertia is directly proportional to 
the cross-sectional area of the links, whereas the link sores— 
ses vary inversely as the thickness of the links. Thus, when 
only Lhu thickness of the links is changed there is no change 
in the maximum link stresses (refer to link stresses in Table 
5.6 for cross-sections of l"x 1" and 1" x 4*', which are 
exactly the same) . However for the same cross-sectional area 
of the links vj^hen the width to thickness ratio is increased 
the link stresses are markedly reduced. This can be obser-^-ed 
by looking at the link stresses given in Table 5.6 for the 
set of cross-sections having the same area, namely, 

1 . 5 " X 1 , 5 " and 2 . 25 " X 1 » or 2 " x 2 " and 4 " x 1 
When the contribution due to KEP inertia towards the input 
end torque is also included the trend is the same as 
diiscussed above as far as the maximum link stresses are 
concerned. The contribution due to the KED inertia towards 
the maximum shaking force, when considered separately, is 
reduced when the width to thickness ratio is increased 
(refer to the maximum shaking force due to KBD inertia axone 
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for the set of cros s-sectious 1.5" x 1,5" and 2.25" x 1 " 
or 2" X 2" and. 4" x 1"). However when the effect of EED 
inertia and HB inertia is considered together there is no 
appreciable difference in the inaxiRiuin shaking force v/ith the 
change in width to the thickness ratio (refer the values of 
the maximum shaking force in fable 5.6 for the set of 
cross-sections 2" x 2 " and 4" x 1"). 

fable 5.7 shows the results of KED analysis for two 
different materials, of which the mechanism links are made, 
fhese two materials are steel and aluminium. Since the 
steel is heavier than aluminium, the maximum shaking force 
and maximum input torq.ue in the case of links made of steel 
are higher than those in the case of links made of aluDiinium. 
fhe maximum stresses in the links made of steel are also highe 
than those in the links made of aluminium. I'his is because 
that the inertia load in the former case is higher than 
that in the latter one. 

fable 5.8 presents a comparison of the results when 
the mechanisms with different link lengths are taken. 

There is a large variation in the link stresses when a 
different mechanism is taken. This table also shows the 
effect of adding a flywheel (having moment of inertia 
= 0.7 Ibs-in-sec") at the output end. It is observed that 
the maximum follower stress is greatly increased by the 
addition of flywheel. In the case of first mechanism there 



























is no effect on the maximum stresses in the crank and the 
coupler due to the addition of flywheel. Whereas for the 
second mechanism the maximum stress in each of the links 
incres-ses with the addition of flywheel. In I’ahle 5»8 the 
analysis is also presented when the input end of the crank 
is instantaneously taken fixed as a cantilever beam. It is 
observed that for the first mechanism the difference bet wen 
the maximum link stresses for the rwo cases, namely, when 
the rigid body degree of freedom is eliminated by the 
method outlined in section 2.7 of Chapter II, and when 
the crank is instantaneously assumed as cantilever, is 
small. However, for the second mechanism the maximum 
stresses in the links when the crank is taken as a canti- 
lever beami are much higher than those when the rigid body 
degree of freedom, is eliminated. 

Figs. 5.4 to 5.6 give an account of the distribution 
of stresses in the crank, coupler and follower, respectively 
versus crank angle for the mechanism A. The damping ratio 
is assumed 0.01 and the stresses are analysed at intervals 
(A0) of 2° up to five cycles of motion of mechanism. Howeve 
the stresses are plotted only for the first two cycles 
except for the coupler stress in Fig. 5*5! which is shown 
for the fifth cycle also. The dotted curves in these 
figures represent the quaSi-static stresses in the respec- 
tive links. It can be observed that the link stresses 




Fig 5.4 Crank Mid-Poini Stres::* vs. Crank Rotation 






Fig 5-6 Follower Mfd-Point Stress vs. Crank Rotation 


determined by complete KEL analysis fluctuate about the 
quasi-suatic stress curve, ihis behaviour is as expected 
in general dynaii/ic response problems and therefore confirms 
the validity of the analysis. The maximum stresses in the 
links obtained by complete IvEI) analysis are greater than 
those obtained by quasi-static analysis. The maximum 
stresses in the former case are more by about 70 r, 114/-' and 
92, -• than those in the latter one for the crank, coupler and 
follower, respectively, during the first cycle of motion, 
i.e. transient state. However, in the steady state the 
difference is only 2;./, 43?® and 8?® between the maximum 


suresses in xlie crank, coupler and, follower, respectively,' 
for the two cases. Thus, except for the coupler stress the 
steady state stresses in the other links, when determined 
using complete ICED analysis, are very close to the quasi- 
static stresses for the mechanism considered. 


Pigs. 5'. 7 to 5.9 represent the distribution of 
stresses in the crank, coupler and follower, respectivelys 
of the same mechanism as above, when the analysis is made 
by aBsuming the crank instantaneously fixed as a cantilevei 
beam, i.e. by considering the mechanism as structure. The 
quasi -static stresses are also shown in these figures by 
dotted curves. A comparative exaitination of these figures 
reveals that the behaviour is better represented when the 
mechanism is analysed by eliminating the rigid body degree 




Fig 57 Crank Mid-Point Stress vs. Crank Rotation (with Crank Input End Fixed ) 









Fig 5.9 Fofiower Mid-Point Stress vs. Crank Rotation (vvilh Crank Input End Fixed) 



of freedom. I'hougli the example cited here is different, the 
result of analysis, when mechanism is considered by eliminat- 
ing its rigid bory degree of freedom, shoM a trend which is 
similar to the one obtained during the experimental deterndna- 
tion of btiesses or strains in the links in some of the 
previous works [34,52]. 

5«2.2 Balancing Analysis and Shaking Porce Kinindzation 

Por mechanism A the shaking force due to rigid body 
inertia only is minimized using counterweights. Iwo circulcr 
weights, with circumference passing through the center of the 
joint, are attached., one to the crank and the other to the 
follower ends, which are supported at the ground pivots. 

I'he radii of these counterweights as well aS angular position 
of their centers are taken as variables. I'he thickness of thes 
counterweights is equal to. the thickness of the respective 
link to which they are attached. Extensions at both the 
ends of the coupler are assumed to the present) the vridth 
of each extension is taken equal to the coupler width but 
the thickness is taken equal to twice the coupler thickness. 
'I'he extension lengths are taken as variables. I'he v?hole 
arrangement is depicted in Pig. 5-10. 

The complete problem is solved as non-linear program- 
ming problem,. In one ca,se the objective function is the 
shaking force. In another case constraints on the maximum 
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shaking moment, the maximum bearing reaction and the maximum 
input end torque are applied and the shaking force is mini- 
mizeo . In a third case the objective function is a combina- 
Cion of the shaking force and the shaking moment, lo make 
d he con bri but ion due to the shaking force and the shaking 
moment tow 9 .rds the objective function approximately equal 
in the beginning, the shaking force is multiplied by a 
weightage factor of 4u . Ihe various results are shown in 
Table 5.9« In this table the results are also given when 
the coupler extensions are not provided and only two circular 
count er weigh! s , one at the crank and the other at the follower 
endSj are used. In both the cases when full force balance is 
achieved the maximum shaking moment, the maximum bearing 
reaction ,and the maximum input torque have all increased 
compared to an unbalanced mechanism. However, by putxing 
constraint on the shaking force, the bearing reaction and 
the input end torque, they can be restricted to be within a 
desired A/'a.lue but at the expense of achieving only partial 
force balance. 

In another example the extensions of the coupler are 
taken ae circular sectors in place of rectangular bars, ihe 
radii and the angles of the sectors are taken as variables. 

Two circular weights are taken a.t the crank and follower enrs 
as before. The center of each circular weight is assumed 

to be at 180° with respect to the line of Joints of the 



















































respective link and the circumference passing through the 
center of the joint (the complete system is shown in Pig. 5.11). 
The radii of these circular weights are taken as variables. 

Ihe op bimizafion resulTS for various cases are given in 
lahle 5 . 10 . Ihe results are similar to ones obtained with 
bar type weights at the coupler ends. 

lable 5.11 shows uhe effect of adding arbitrary 
masses at the ends of the links of the mechanism A on the 
stresses in the links and on the contribution of IvBD inertia 
towards the shaking force etc. The stress analysis is 
made with an assumed value of damping ratio of 0.05 and 
angular interval is taken as 5 ^ > The first row contaans the 
maximum link stresses and the maximum shaking force etc. 
when there are no count erweighcs attached to the ends of the 
links. The rows that follow the first row show the maximuir.. 
stresses and the maximum shaking force etc. when a counter- 
weight is attached in turn to the end of each link. The 
counterweight attached to the coupler ends is circular 
sector having radius as indicated in the table and the 
angle of the sector is one radian. The counterweight plactu 
at the crank or the follower end. is circular with its 
center at 180° with respect to the line of joints and 
circumference passing through the center of the joint. Tne 
addition of counterweight at the crank end produces a reduc- 
tion in the shaking' force, particularly that one which is 
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due to rigid body (fiB) inertia, but the maximum stress in 
the crank is increased. The addition of counterweight at the 
follower end causes a reduction in the maximum shaking force 
due to rigid body inertia but the maximum shaking force, 
the maxinium shaking moment and the msxiitium input torc^ue due 
to llED inertia are increased. The maximum stress in the 
follower is also increased. The addition of counterweights 
of the amount indicated at the coupler ends has adverse 
effect on balancing of the mechanism hut it produces a 
reduction in the maximum stress in the coupler. 

5.3 Optimization Results 
5 . 3.1 First Problem 

The example, chosen for this problem is a crank 
rocker mechanism to guide a rigid body through three given 
positions. The following data is assumed; 

V/hen the crank is turned through 30° and 90° from 
some initial reference position, the rigid body is displaced. 
in the x~direction and y-direction by an amount -7.0" and 
2.0" and -20.0" and -4.0", respectively, from some initial 
position. The initial cross-section of the crank is taken 
as 3" X 1" and that of the coupler and follower as 2" x 1". 
The cross-sectional area of the each link is varied by 
taking a variable in the design vector which is used to- 
divid.e the ini-cial width of the link to determine width at 
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a new design point. Ihe width to the thickness ratio is 
maiiit; e.inecl as 2 unless the thickness is reduced to the 
extent of 0.2", when it is kept constant at 0.2" and in this 
case only bhe width is allowed to vary. Ihe upper bound on 
c he iBaximun link stresses is taken as 1500 Iby^sfi in. Oince 
the coupler point, where the body to be guided is attached, 
comes out to be very close to the line of joints of the 
coupler (distance Ig in kig. 5.12) initially as well in the 
final design for this example, the effect of coupler extension 
is neglected during the liED analysis to reduce the compuxa- 
tional time. 

Table 5.12 shows the initial design and the optinmm 
design with and without including the effect of kEL inertia. ■ 
The objective function in both the cases is the shaking 
force. The maximum shaking force, the maximum shaking moment 
and the maximum input end torque are reduced by approximately 
l/15th, l/8th and 1/llth times the original values, res- 
pectively, for the case in which the contribution due to 
KED inertia is Included along with the contribution due to 
rigid body inertia towards the shaking force etc. This 
reduction is- partly achieved due to a change in the link 
lengths and partly due to a reduction in the link cross- 
sectional areas. The stress, constraints are far from being 
active. Here apparently, one may think that since thc- 
link stresses are small compared to the allowable limit, 




Table 5.12 Optimization Resuits-Problern uajancKb) 
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on© may reduce the area of cross-section further and achieve 
a further reduction in the shaking force, ho doubt reducing 
the lin..h cioss— sectional areas further may reduce the 
contiibution cue to the rigid, body inertia towards the shak- 
ing 1 01 ce etc., but the contribution due to KED inertia 
towards the shaking force may increase due to high elastic 
accelerations because of slender links, fhis fact is 
depicted in I'able 5*13 where results are given when the 
link cross-sectional areas are changed arbitrarily at the 
optimum design point, i'here appears to be two local tiniUo^ 
close to each other. I'he significance of these local miniiua 
can be explained as follows. Around the first local minimum 
the reduction in the shaking force due to rigid body (SB) 
inertia is less compared to the increase in it due to the 
KSL inertia when link cross-sectional areas are reduced 
at first. However, the reduction in the shaking force due 
to RB inertia again over compensates the increase in the 
shaking force due to iiED inertia when link cross-sectional 
areas are further reduced about the first local minimum 
giving another minimum close to the first one. Around the 
second local minimum the increase in the shaking force due 
to 'EllD inertia, when link cross-sectional areas are reduced, 
is more pronounced and the net shaking force increases for 
ever. It can also be observed from fable 5.15 bhat the 
maximum link stresses increas.e with decrease in the link 



Table 5.13 Effect of Arbitrary Variation in Links Areas at Final Optimum Design 

Problem l(a ) 



Q 

^ S IH ^ 

i oP 

< P O 

CL t~ 


{£) cn 

CD 

CO O 


CM 



. ’O 

cn 

. O' 

m 

00 

. 0 


i>4 

•'■1 



r^' 

r-j 

01 

r'.i 


P'v' 


uo 

T ' 

cn 

cO : : 

.C"" 

f- 

0 

■ IX’J 1 

t>F 


cn 

CO 


■N 

m 

n- 

CD 

cn 

cn 

CO 


(T-' 

fO 

CD ■ 

t — 

CM : 

nJ 

CM 

CM 

CM 

CV4 

(SI 

<NJ 

r-4 

rr^ 

rn 

■•D 

-t 

■vD ■ 

■ ' -4 ■ 

vT 

■,T ■ 

“C- ■ 

■ -^7 ■ 

. --T 



00 00 

CO CjD 


i 

■ DO ' ■: 

1 

ro 

1 

00 

. ^ i 

fO 

1 

i 

■ -r 

1 

■t 

1 

-t 

1 

-T 

i 

tr^j ■ 

1 

0 

C'l 

CD 

CO 

in 

CD 

■X 


■ <0 

• 0 

■ r^: 

CO 

CO 

00 

CD 

4* - 

CD 

in 

r— 

CD 

T". 

C" 

CO 

in 

in 

■J- 

•4- 

ro 

CO 

CO 

CM 

r— 

c: 

r— 


-r- 

r— 

'■■V 1 

r- 

r~ 



r“ 


»—• ■ 


0... 

”> 2- Z 

^ f LU c 

X <i 2 ■ ”• r^4 

< X O ,a 
2 : 1 /) iS'-' 


un vj- 

,n . ^ 


, [v. <0 -sj CO CO o c 4 a c 

00 CM cn cD cD ^ ^ 

on cNj r- CD <0 cn iT) cn r-, r- 

(V) ro 00 cNi 04 O'! ri n -o •T 


SO , 

”3 2 UJ D' 

2 X O [H., 

” cc XI 
20 ^ 


r-H CO 

<0 -j' 

on ; fO 


<0 in 
i . cO 


'cn Vo ■ O' , ^ ''' 

^ CO ^ F* CO 


' ,Z. ■ 

0 ce cNt 

< P o OP 

UJ (J 4-~ . 

q: o y o 

<c '■Q,: < "■ 

LU 

'.OC..', 


(D f o . ' OV;. 0 , C3 


O ■ 


m CO 


o 0,0 o 


o o 

00 CD 


^ c» CD 

000 


CO, : 'o 

o cn 



^ Q ' S S 0 ^ 

< £ a t5 P X 

fc ^ ^ ^ o 5 

I s 3 s fi ^ 

5. > '” cc 


ax - 5 : 

u S S ? 

^§3 2 

2 ^ 2 : qI 

^ ^ o 








cross-sectional areas. 

P grabs oi optimization process 
rs s.o.n in lalle 5 . 14 . - can la cos ervea t. at r.e process 

is first converged to an optimnn. design .irlout an, signi- 
frcant change in the linh cross-sectional areas. One of the 
Grashox constraints becomes active, xhus at first the reduc- 
tion in the Shaking force is only due to a change in the 
link lengths. At this design point arbitrary variation in 
the area reduction factor ie made and the effect is 
observed. I'he results, which are shown in fable 5.15, 
indicate that a further reduction in shaking force is 
possible by a reduction in the link cross-section areas. 
Therefore, the problem is further optimized starting from the 
first optimum design, point. However to save computational 
time the link cross-sectional areas at the initial point 
are arbitrarily reduced to halves of their values at the 
frrst optimum design point. The progress of this optimization 
process is also shown in fable 5.14. I'lie fact that optimum 
design thus obtained is a final one has already been explained 
with the help of results shown in fable 5.13. 


When the effect of KBD inertia is not taken into 
account the maximum shaking force, the maximum shaking 
moment and the maximum input torque continue to decrease 
till, one or more stress constraints become active. The 
progress of optimization process is shown in fable 5.16. 
Here also it can be observed that at first the reduction 
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Table 5-16 Details of Optimization Progress - Problem I (bj 


1 


; or 







riri.| 

UL) 








: 5 


0 

0 

K 

CNJ 

r— ■ 

CO 

o 

’ J 

(N 

to 

01 

0 

0 

,'Mg 

00 


cn 

; 

U) 





!./) O 




CM 




Ll 

t/).E 




CM 

'-I" 

n 

c/) a ry 

f; 









cn 

.CM 


vf'. 

£■ 

in 

if) 


: ■■ 

'/CO 

8 2 

a. 

a-: .. 

<jD 

0 . 

■ O 

0^ 

\o 

o 

D2 : U 

2 : 1 .-.-... 

T;. 

• ' ’ ' ■' 

r 

r- 

■ - J ■ 

... 

cn 

1 

>'0 

CD 

. . 

( 

^ -y 

CO 


CM 

Mf 

CM 

<r 

0 

2 Z 

' 


1 

cn 

M- 

cn 

r" 

CN 

Q: 

O’ 

CD 

iT) 

' o> 

■f-' 

r*' 

CM 

'T 

CD 

■I 

CO 

■t 




i , 





r~ 

. 

CN 

! 

■ 

0 

0 

0 ■ 

cn 

or 

• t— 


0 

CO 

MT 

n 


in 


CO 

CD 

vT 

CM 

a) 

mcc 

6 

0 

0 

0 


r—~ 

r*". 

■ 

p-‘ [ 1 V_^, 


0 

0 

0 


cn 

CD 


cn 

1 CO 

CO 

CO 

CM 

fO 

cX) 

< X ‘i) 


in 

tn 

M- 

m 

n 

sT 

ct “ 

0 

. 0 ' 

6 

0 

0 

0 

0 

t- < • 







^/)C!C . 

CO 


CD 

cn 

cn 



z r: 

OD 

CM 

r- 

, CM 


r— 


O2 

•<r 

fNi 

CM 

<0 . 

m 

n 

nJ 

cj 5 

6 

0 

6 

6 

0 

0 

0 









^ w t 

— § ^ 0 

r~ 

0 

CM 

CO- 

CM 

CO 

r-H 

; 

00 

.(■X 

CO 

CM 

- 4 - 

H- 2 < 

{N 

n 

CO 

CO 

0 

. 

.'O 

Z 

0 

r— 

iD 

a 

0 

in 

n 

'■ C'.V 

^ 5 

C^J 


■ : — 

fSl 

'CM 


"'.J' 

-cr 

<r 

m 

0 

CO 

r-“'. 


3 

CM 

' r— 






u. 








^ g 

cn 

■■ 

0 

CM 

: r-T ' ' 

CO 

0 

r— • ■ 

0 t- 

CD 


0 

■■■ »•■ ' 

CM 

n 

fO' : 

UJ 0 

<N 

. fO 

0 . 

CJ 











m S 

0 LL 

<— 



"■ ■ ' 

|_ 



AREA 

EDUCT- 

ION 

FACTOR 

5 .00 

5.1 6 

0 

0 

V— 

1,94 

n 

04 

r* 

cn 

CO 

■ XT' 

0 

s>'.. 

QC 








^ >: ^ ' 

f- L d; 

""H 

'-■o- 

' , ' ' ' 

■■■. 0 . 

o- 

' 0^ 

0 

■\,0 :: 

■ ■' 

■d S liJ 
< a H 

0 

0 

0 


■■■O' 

■ 'n ■■ 

" 0,2' 

Z < IL) 

tn 

in 

m 

LD 

iT) 



. Gt - 














■H-- 



< 

■ i' ■ 



■ '•• , 

Q: ■ 

O 

Q'''' 

0 : 

x:';.. 

TIMIZ 

iON 

STEP 

< 

p 

z 

.'h^ ■ 

tn 

S®:: 

< 

■ it— ■ 

2/ 

FIRST 

< 

L- 

if) 

LU 

c 

C) 

U! 

a. 


U- 

*-«* • ■ 

CC 

to 


0 











12G 


in shaking.' force is dup tn 

. bn ti,e link lengths and 

there xs no appreclehle =hen«e In the llnh oroes-eecthonel 

areas. Phe area redaction factor le then scaled dc» and the 

prohlet, IS further optimised. ,t the final „ptl„t. design 

point the llnh stresses are close to the allowable ll„it . 

Obsem/ing the trend from the tablp it n v 

able lb Can be expected that 

any farther reduction in areas of the linhs ,nll cause the 

llnlc stresses to cross the allowable licit, fhe link areas 

at the optimum design, point are much smaller in this case 

than those in which the effect due to Ki'D inertia is also 
included . 


5 •3.2 Second Problem, 

I'ho shaking force is minimized for uiechenism A when 
it is running ac 300 rpm, The initial cross-section of the 
crank is ta.ken as 1.5 "x 0.5" and that of the coupler and 
the follower as l.O" x 0.5". two circular counterweights 
are ta.ken one at the crank and the other at the folloiver end 
with their centres at 180° v/ith respecr to the line of pivots 
of the respective links. The circumference of the counter- 
weights passes through the center of the respective fixed 
pivots. iwo Sector type counterweights are taken one each at 
the two ends of the coupler. The complete arrangement is 
indicated, in i'ig. 5.11. Ihe design variables are the radii 
of t,h8 four weights and the angles of the sectors. The 
variation in the link areas is achieved by, taking a variable 



in the design v3cxor dividing the width of the each link 
by tiiis variable at each design point. ihe ratio of the 
width to the thiclmess is maintained at 2 till the thickness 
IB reduced to 0.2" at which stage it is taken as a constant 
(0.2" ) and only the width is allowed to vary, ihe thickness 
of the each counterweight is taken eiual to the thicicness of 
the link bo which it is attached, hon-negativity constraints 
on the sector angles are imposed. Also the sector angles 
are not allowed to increase beyond 271; by putting constraints, 
ills upper bound on the maximum stresses in the links is taken 
as 1500 psi. I'he initial radii of the counterweights are 
taken as 2" and the sector angles as one radian. 

I'lie results of optimization are shown in iable 5 "17 
with and without KSD inertia effects. I’he objective function 
in both the cases is the maxiro.uffi shaking force. The reduc- 
tion achieved at the optimum design in the maximuin shaking 
force, the maximum shaking moment and the maximum input 
torque is approximately 80, 5b/'’ and 667^ of their original 
values, respe ctiAT'ely, when the contribution due to nEI) 
inertia is taken into account. The change in the radii of 
the counterweights is very little (except the radius of the 
counterweight attached to the, crank end) at the optimum 
design comipared to their values at the initial design point. 
However, the reduction in shaking force at the optimum design 
point Cannot be claimed only due to the reduction in cross 
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sectional areas of the links, fhis fact is made clear by 
determining the maximum shaking force, the maximum shaking 
moment ano. tne ma.xim'um stresses in the links when the link 
cross -be ctiona.1 areas are reduced arbitrarily without any 
counterweight attached to the links. I'he results are shown 
in iabl© 5 « 18 . It can be observed that the shaking force does 
not decrease below a value of about 14u lb. However, with 
counterweights the optimum value of the shaking force attained 
is 56 lb. Any further reduction in the link cross-sectional 
areas a,t the optimim design point does not produce a reduction 
in the maximum' shaking force. In fact the shaking force ever 
increases when the link cross-sectional areas are reduced from 
the optimum value with the optimum set of counterweights. 

I'hese results are shown in I'able 5.19. I'he progress of 
optimization progress for this case is shown in 'I'ahle 5.20. 

In the Second case, in which the contribution due to 
KED inertia, towards the shaking force is not included, the 
reduction in the maximum shaking force is achieved partly 
due to the adjustment in the counterweight parameters and 
partly due to a reduction in the link areas. The progress 
of optimization process is shown in Table 5.20. It can be ooser 
ved that the link stresses are approaching to their constraint 
boundries near the optimum, design point. Thus a furthex 

A 

reduction in the link areas may cause a reduction in the 
shaking force but the link stresses may cross the upper 



Tabfe S-TB Effect of Arbitrary Variation in Links Areas of Original Mechanism 
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Table 5-20 Details of Optimization Progress - Problem 11(b) 










limit. The coun..er„-eieht at the crank end reduces to nearly 
aero, which is contrary to the earlier obserrations shown in 
'X'ab3.es 5.9 anc 5.10, where the shaking force due to rigid 
body inertia only is niinimized without stress constraints. 

Inis is because tiia,t Ihe counterweight at the crank end 
causes tin increase in iis stressj on which a constraint' 
has been imiiosed, thus it does not allow a counterweight 
to be present ac- that end. As such, it is not possible to 
achieve fiill force balance bj. redistribution of link masses 
when stres;.; constraints are present. 

In tlie la.st ro'w of fable 5.17, the results are 
presented when the counterweights, determined according to 
the method presented by Berkof and Lowen [b], are attached 
to the crank end the coupler ends to achieve full force balance 
I'he meximarn shaking iorce becomes practically zero but the 
maximum shaking moiuent and maximum input torque are 20‘;i and 
40/i more than thoir values at the optimum design point 
obtained using the present procedure. Moreover the maximum 
stress in the .follower it very high and has surpassed the 


allowable limit . 



CHAPl'ER VI 


OOiMULUblONo Afflj PROPOSALS POR PTJi'URL WORK 


2.1 Conclusions 


In the present work only one example in each of 'uhe 
two types of the pro Declares to achieve balance in a planar 
mechanisin has been solved. k>ince it takes a long compaua" 
tiona.L tinie to solve each problem (about 6u min- of CPU 
time on the average on DIG 1090), it is not feasiole to 
solve many examples. However, the results of the numerical 
example in tli::; each category can be interpreted to give 
general conclusions for not only a four bar planar mechanism 
but for any planar mechanism in general. ibe fol ' ' B 
conclusions are drawn from the results of the, two numerical 


examples solved . 

(1) It is extremely important to taie into account tns 
contribution clue to KED inertia towards tbe shaking 
force anc. the shaking moment to achieve a balance in 


mechanisrcs with elastic links in a real sense. 

2) Ihe inclusion of the effect of the inertia, limrue 

the amount of balance that can be achieved. 

• t oHTieve full force or/and full 

3) It is nox possiDle to achie 

• V 1 -r-inn rtf the link masses 
mom. 3 nt balance by redistri 
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because the contribution due to KBD inertia towards 
the shskiiin iorce and the shaking moment cannot be 
completely balanced by any combination of counterweights 
However, a reduction in the shaking force and/or 
shaking moment due to KBD inertia is possible by using 
s. suita.ble set of the counterweights. 

(4) When the effect of kBD inertia is not included and the 
me cliani, Giv: it. lully force balanced using suitable 
counterweights the dynamic sti-esses produced in the 
links are xn;r,y high because of the increased inertia 

01. the J.ink.ti, IMiut! a kin etoelasto dynamic analysis must 
be made to deteCTlne the link stresses when a mechanism 
with ela.Btic links is to be fully force balanced using 
counterweights even if the contribution due to lu^l 
inertia towards the shaking force is neglected. 

(5) It was 0 St, a blit, hod by Conte and his group [15] that trie 
shaking iurce, shaking moment and input end. torque d-.e 
to rigid body inertia forces' could be reduced by 
al'bering, tliC; kinematic design. I'he present stucy 
reveals that this procedure is also effective in 
reducing the contribution due to KBD inertia towaroo cni: 
.shaking force, shaking moment and input end torquu- 

(6) 'i'he redaction in the link areas is very effective, ii. 
minimising the unbalance in the mechanisms . However, 
the liiik stresses must ue deterid.ned at eacix stabW 
check that they do not cross' the allowable limi^ • 



( 7 ) 


When the contribution due to KEh inertia towards the 
shaking force and the shaking moment is not included, 
the link cross-sectional areas continue to decrease 
with the unbalance in the mechanism becoming smalj-er 
and smaller till one or more of the link stresses 


( 8 ) 


( 9 ) 


approach the allowable limit, 

The assumption made in determining the gradien'Cs oi o 
various function as outlined in Chapter t has vxorkec- 
well . It has reduced the computational time considei 
ably for the second problem in which the maximum 
values of the various functions occur mostly curing t 
first quarter part of the motion cycle Oj. the mecxiaiii 
However, for the first problem, at the latter stages 
optimization process, the maximum values of some of 
the functions occur during the last quarter of the 
motion cycle of the mechanism therefore the analysis 
has to be carried out till the end of the complete 


ycle of motion for most of the time- thus there is 

ery little saving in computational time- in this casv. 

‘he use of the gradients of the link stresses and 

maximum shaking force to predict them at a new desi.r 

,oint during one dimensional minimization has success 

1 A also resulted in a consioe.'’, 

iully worked, inis n«s aisu x j 

saving in the computational time. 
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2.2.2 becond Category 

Under the second category the following suggestions 
are proposed. 

(a) lit the present work only one design variable is used 
to make a change in the area of all the links to save 
the computational time by limiting the number of design 
Trariables as far, as possible. However, the maximum 
stress developed in a link and. the contribution due to 


KEt) inertia of a link towards the shaking force and 


oha.king, li.oment are different for different links, 
iherefore to be more realistic, the variation in the 
area of each link must be independently made. Moreover, 
the ratio of the width to the thickness of a link is 


taken a constant (2 in the pres-nt study) . ' For a 
mechanism of given link lengths and input crank speed 
the rigid body inertia forces depend only on the link 
cross-sectional areas and. are independent of the shape 
of the link cross-sections. On the other hand, ^he 
KIIID inertia forces depend on the shape of the link 
cross-sections since the elastic deformations and as 
such the elastic accelerations depend on them. Thus, 
the width to the thickness ratio of each link ma^r also 
te allowed to vary or soB.e other shape of the link 
croes-sec-clons may be chosen to arrive at a better ^ 


design . 



2.2.2 Second Category- 

Under the second category the following suggestions 

are proposed. 

(a) In. the present work only one design variable is used 
to make a change in the area of all the links to save 
the computational time by limiting the number of design 
variables as far. as possible. However, the maximum 
stress developed in a link and the contribution due to 
KBD inertia of a link towards the shaking force and 
shaking moment are different for different links. 
I'herefore to be more realistic, the variation in the 
area of each link must be independently made. Moreover, 
the ratio of the width to the thickness of a link is 
taken a constant (2 in the pres':-nt study). ' For a 
mechanism of given link lengths and input crank speed 
the rigid, body inertia forces depend only on the link 
cross-sectional areas and. are independent of the shape 
of the link cross-sections. On the other hand, the 
KED inertia forces depend on the shape of the link 
cross-sections since the ela^stic d.ef ormations and as 
such the elastic accelerations depend on them, ihus, 
the width to the thickness ratio of each link may also 
be allowed to vary or some other shape of the link 
cross-sections may be chosen to anrive at a better 
design. 
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(b) lh.e KBD ajialysifci is made by considering the contribution 
due only to rigid body inertia forces towards the load 
vector, i'he additional inertia forces produced due to 
vibration of the links are not included in the loao 
vector for displacement analysis. A more realistic 
approach be made by taking into account the contributio 
due to vibrational inertia forces towards the load 
vector used in making KEl analysis. 

2,2.3 Third Category 

ihe effect of ICE® inertia forces on the unbalance m 
other planar mechanisms may be imz-estigat ed and. the 
tiveness of the present procedure in arriving at at P ' 
design be d.eterniined . 
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APPEEDIX A 

eemema'I'ic analysis 


A kinematic analysis is made to obtain positions^ 
angular velocities and accelerations of the various links 
in a four bar mechanism when their values for input crank 
are given. 


Referring to Eig. D.l(a), let 1^, 1 ^, 1 ^ and 1 ^ are 
the lengths of the crank, the coupler, the follov/er and 
the ground link, respectively, measured between their joint 
centers, 6^, 0^ and are the angular positions of the 
crank, the coupler and the follower, respectively, measured 
in the counterclockwise direction with respect to X-axis 
as shown in Pig. D.l(a). The angular positions of the 
coupler and the follower are related to the crank position 
and link lengths as follows [ 53 ], 


CO s 0 2 = ( A P + C j Cg C sin 6 /D 

sin 0 2 = (A C^C-CjB sin e^)/D 

cos = -B cos 02 + C sin 

sin 0^ = -B sin $2 - ^ 6 ^ 

where A = 1 - cos 0 ^^ 

= C4-C5 cos 

= ^[1 - ( '^1+^2 ^ 1^^^ 


( A-1) 


C 


{A- 2 ) 

. . , Contd. 



IdO 


and 


D = - 2C3 cos 0^ +1 

B = + Gg cos 0 ^ 

0l = (^2^ + - l4")/2l2l3 

^4 “ ’*’^2 ~^3 '''^4 ^ /^^ 2^4 

^5 ^1/^2 

^6 = l 3/'-4 


( A- 2 ) 


( A- 3 ) 


When the input crank runs at const axit angular speed the 
angular velocities and accelerations of the coupler and the 
follower are given by 

aj_ = w- 1, sin ei-0x)AoSin (B -Q r) 

2 J- ^ (A-A) 

= “1 ^1 (0 2_“0 2 ) A 5 ® i “- 

“2 " ~*'^l“l^ (01-03) + i2“2^ (02-^3) l 3 “ 5 ^)/ 

“1^ sin (®2~®3) 

2 

a-^' - cos (®i“0 2) ^3*^3^ (62-03) ^ 

1^ sin (8 2 -0^) (A-' 5 ) 

where and are the angular velocities of the crank, 

the coupler and the follower, respectively and a2 and are 
the angular accelerations of the. coupler and the follower, 
respectively. 



IPPEKDIX B 


element ma't rices 


'The element stiffness matrix k, mass matrix m and the 
load vector p, which are referred in Chapter II (eqs. (,2.21) 
to (2.23)), for an element of uniform cross-section in element 
oriented coordinates system are given by 


k = 


1 

l^(l+9) 



I 


1 0 

1 

12BI 



i “ 

bBIl 

(4+9)Bir^ 

SymiT'etric 

~EA1^ 0 

0 BAl^ 


0 

-12EI 

-6EI1 0 

12EI 

0 

6BI1 

(2-9)BI1^ C 

-6EI1 (4+9)BI1 


m = 


P A1 

840 ( 1+cp ) 


f 280 ( 1+9 )^ 

0 ( 312+588q)+280(p^) 

0 (44+779+359^)1 

140(1+9)^ 0 

0 (108+2529+1409^) 

0 -(26+639+359^)1 


(B-1) 


SymiTietri c 


(8+1409+ 79 '"')!^ 

0 

( 26+ 659 + 359^)1 
-( 6 + 149 + 79 ^) 1 ^ 


Symmetric 


280(1+9)^ 

0 ' (312+5889+2809^) 

0 -(44+779+359^)1 (8+1409+79^)1^ 

(B-2) 
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P AL 

P - " 120(1+^ 


I 20(l+cp)(2a^+ag^) 

( 42+40cp ) 3-j^y +( 18+20q)) 
(6+5cp)l + ( 4 + 59)1 

20( 1+cp) ( s..^x+2a.j^^) 

(18+209) a^y + (42+409)agy 
-(4+59)1 + (6+59)1 a^y 


(B-3) 


where 1, A and 1 are the length, cross-sectional area and 
second moment of the cross-sectional area, respectively, of 
the element (link) . p and E are the density and modulus of 
elasticity of the element material, a^y, a^^^ and a^y 

are components of the accelerations in the x and y directions 
of the two ends A and B of the element, respectively, as 
indicated in Eig. 2.2(c) of Chapter II. 



APPBffl)IX C 


ABSTRACT OR SOlUl'ION PROCEDURE ROR BIGEIWALUE PROBLEM 


To determine the eigenvalues and eigenvectors of the 
dynamical matrix of the eigenvalue problem ( eq. (2.47)) given 
in Chapter II an algorithm due to G-rad and Brebner [45] has 
been used.. This algorithm gives all the eigenvalues and 
eigenvectors of a real general matrix. It involves QR double- 
step method to compute the eigenvalues and inverse iteration 
process is used to find out the corresponding eigenvectors. 

A brief summary of the procedure adopted in writing the 
algoritlmi is given below. Ror other details of the algoritlim 
the original paper mentioned abo%’e can be referred. 

To obtain a better accuracy in the results the follow- 
ing two modifications are. first made. (i) The absolute sums 
of corresponding rows and columns are made approximately 
equal by scaling the matrix using a sequence of similarity , 
transformations. (ii) The resulting matrix is normalized, 
so that the value of Euclidean norm is equal to one. 

'The matrix is reduced, to an uppsr-Hessonberg form using 
similarity transformations (Householder’s method). Using QR 
double-step iterative process the He-ssenberg matrix is reduced 
to a form so that all elements of the subdiagonal that converg 
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■to zero are in modulus less than 2 |lH|lg, where t is the 
number of significant digits in the mantissa of binary 
floating-point number. I'he eigenvalues are then obtained 
from the reduced form. 

Inverse iteration is applied to the upper-Hessenberg 
matrix until the absolute value of the greatest component 
of the right hand side vector is greater than the bound 
2'^/(100N)5 where N is the order of matrix. Generally an 
additional step is performed to get the computed eigenvector 
after achieving the above bound. However, the residuals are 
determined at each step and if at any step they are greater 
in absolute value than that of the previous step, then the 
vector of the previous step is taken as the computer eigen- 
vector . 



APPENDIX D 


EXPLICIT EQUAITOKS POfi DYNAMO JORGE ANlEYblB 


D.l Introc'.uctioii 

I'here a.re many methods availsble for the dynamic 
force analysis. Out of which most commonly used methods 
are directly based on Nevrton's second law or principle of 
Arirtual work [55,56 ]. hatrix methods [57,58] are also now 
available for this purpose. However, since in the present 
study a simple four-bar mechanism, has been chosen for 
numerical examples, the dynamic force analysis is made by 
applying Newton’s second law to the free body diagrams of 
the links, which readily gives a set of equations. The 
shaking force, the shaking moment, the input end torque and 
the bearing reactions are then determined using these equa- 
tions. The complete procedure is presented in section D.2. 
Section D.3 deals with the procedure to determine equivalent 
dynaDlcal system when counterweights are attached to the 
links . 

D.2 Dynamic Force Analysis 

In Jig. D.l (a) a planar four bar mechanism at an 
instant of time is shown. The links 1, 2, 3 and 4 ale 
jointed together of which link 4 is fixed. when the m.echaniBm 




is in motion the links are subjected to angular and linear 
accelerations which vary with respect to the mechanism 
configuration. 'I'his gives rise to varying rigid body 
inertia forces v^hich in turn produce varying reactions at 
the ground bearings (supports). 

To determine these reactions at any instant of time 
we consider the free body diagrams of the three movable links 
1, 2 and 3 in the positions corresponding to the instant 
considered as shown in Pig. D.l(b), 

Let 1 ^ be the- distance between tvK) joints of ith link, 

6 ^ be the angular position of the link, r^ the distance of 
the centroid of the link from the left end joint for the crank 
and the coupler whereas from the right end joint for the 
follower, be the angular position of the centroid with 
respect to the line of joints of the link. The angles are 
measured positive in the anticlockwise sense. R- ■ and R^ , 
denote the X and Y components of the forces exerted by the 
jth link on the ith link at the joint. DXi, DYi and are 
the components of inertia forces in the X and the Y direc- 
tions and inertial torque of the ith link, respectively. 

Tj^ is the input torque at the driving crank. Then from the 
free body diagrams of the links (Pig. D.l(b)) we have the 
following set of equilibrium equations. 



1 


iS 



+ 

■^12 

+ 

DXl 

= U 


+ 


+ 

DYl 

= 0 

^21 

+ 


+ 

DX2 

= 0 

^21 

+ 

Ry 

^23 

+ 

DY2 

= 0 

^32 

+ 

li? 

^bi 

+ 

DX 3 

= G 

^32 

+ 

Ry 

R34 

+ 

LY 3 

= 0 


'^'l + *^1 “ ^12 ^12 ^1 “ (D“l) 

__ DXl-r^ Bin (® ■*■ 9 i) + 011 cos = '^' 

^2 ” ^23 ^2 ®2 ^23 ^2 ®2 

- 7 JX^-r^ sin (Q 2+^2^ DT2*r2Cos (0 2’*'<P2^ “ ^ 

C3 “ 1^ sin I2 cos 0^ 

~ DX3'BG-sin (Q^->p^) + BY3’ BC- cos(0 = U 


I'he new parameters BC and used in the last of the eq.s. 
(D- 1 ) are defined as follov;s from triangle ABC shown in the 
third, link [Pig. D.l(ti)]. 


and 


BG^ = 1^ + r| - 21^r2 cos (p^ 

-1 ^3 

'i’ ^ = sin sin cp^) . 


(D-2) 


Keeping in mind the following identities the eqs. 

(D-1) relates the nine unknown quantities, i.e. the reactions 
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at the four joints in the two directions and the input 
torque. 


X 

^21 

%2 

- -B 2 I 

^23 


^^23 

- -R 32 

^34 

- -^43 


- -^^43 

»41 

pX 

- “^14 


^14 


(D-3) 


ihus using identities (h-q) in eqs . (E-1) one gets nine 
linear non~homogenous equations in nine unknomis which are 
solved to give these unknowns. 

If HU be the mass of the ith link, I^ the moment of 
inertia about center of gravity of the ith link, 
and are the components of the linear accelerations of the 
centroid in the X and Y directions and the angular accelera- 
tion, respectively, of the ith link, then one has 


L'Xi = 

-ki 

axV^ 

DYi = 

-m. 

ap 

Ci = 

-li 

“1 


(D-4) 
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'i^xe positions of the links, their angular ''elocities and 
accelerations for a given position and speed of the input 
crank is obtained from the kinematic analysis of the mechanism' 
presented in Appendix A. the linear accelerations and 

aYj^ of the links (i = 1,2,3) are then given by 


aX^ 

= 

C\J 

1 — 1 

3 

1 — t 

1 

cos (0^+(p^) + sin (0 2+cp2^)] 


aYj_ 

= 


sin (0 ^+.:p^) - cos (®j_+rp2_)] 


aXj 


-rgtiDg 

P 

cos(^ 2+92^ ''’“2 

01 

afg 

= 

-rpo./ 

sin( 0 2+92^ "“2 2'^'^2^ ^ ”^1^1 

01 


= 

2 

l-r.^r 

O'- 3 

cos(6^+cp^) + a.-^ sin (6^+9^)] 


aTj 

= 


sin(0 ^+9^) - cos (0 ■5+9-^)] 

(D-5) 

where 0 

and 

are 

the angular velocity and the angular 



acceleration of the ith moving link. 

■file components and of the shaking force in the 
X and Y directions , respectively, are given by 


- ®^4i ^43 




- ^^41 + ^43 


(D-6) 


Using first six of the eq.s. (D-l) , identities (D~3) and eqs. 
(D-6) one gets 


FS^ 

= DXl + DX2 + I>Z3 


FS^ 

= DYl + DY2 + DY3 

(D-7) 

FS 

= V7_(FS^)^ + (FS^)^] 
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where FiJ denotes the net shaking force acting on the founda- 
tions . 

If the X-axis is taken parallel to the ground link 
(i.e, fixed link 4) the shaking moment about a point midway 
between the two fixed supports due to ground bearing reac- 
tions is given by 

= 1 (B-8) 

where, 1 q is the distance between the fixed supports. The ■ 
net shaking moment transmitted to the foundations is given by 

m = MS' - Tj_ (D- 9 ) 

The bearing reactions at the four joints are given by 


Bill = 

BR2 = + (^ 12 )^] 

BB 3 = 

BR 4 = V[(R34)^ + (^34)^] 


(D-10) 


The maximum values of the shaking force FS, shaking moment 
SM, input end torque Tj_ and the bearing reaction during the 
entire cycle of motion of the mechanism is determined by 
performing the above analysis for various configurations of 
the mechanism differing by a preselected angular interval 
(A6 ) of the input crank from 0° to 360 • 
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D.3 Determination of Equivalent Dynanucal System 

Since a four bar mechanism is chosen for all the 
numerical problems solved in the present study, the deter- 
mination of equivalent dynamical system is presented only 
with reference to the four bar mechanism. 


Consider a mechanism with counterweights attached to 
its links as shown in Pig. 5*10 of Chapter V. ihe complete 
elaboration of the arrangement is also given in the same 
chapter. Let and are the radii of the circular counter- 

weights attached one to the crank and the other to the 
follower, respectively a and Up are the angular positions 
of the counterweights with respect to the line of joints of 
the respective links. I'he coupler extension lengths are Ij^ 
and 1^ corresponding to the left and right ends of the 
coupler, respectively. Let and are the mass and the 
moment of inertia of the ith link about its centroid without 
counterweights, and m^ and 1^^ are the mass and the moment of 
inertia of the ith link about its centroid with count er- 
wsight( s) , respectively. 


The mass and the moment of inertia of the each circular 
counterweight about its center are given by 


m 


A 


m-. 


nR/ 



B 


P 

P 


(D-ll) 



(D-12) 


T 1 2 

Ia - 2 \ 


where t , m and I denote the thickness, mass and moment of 
inertia of the count erweighu and the suffix A and B stand 
for the crank and the coupler counterweights, respectively, 
p is the density of the material of the counterweights. 


The position of the centroid of the link and. its 
attached counterweight combined, in a coordinate system in 
which x~axis is taken along the line of joints of the 
respective link and y-axis perpendicular to it and the origin 
is . the center of the left end joint for the crank and the 
coupler, and that of the right end joint for the follower, 
is given by 


Xj_ = cos aj^)/(in’+m^) 

y^ = ( sin a^) /( m' +m^) 


(D-13) 


and 

Xj = (Im^l ^ + mgRg cos ag)/{m'+mg) 

(D-14) 

y^ = sin ag)/(m' + m^) 

where x^, y^ and x^, y^ are the coordinates of the centroids 
of the crank and the coupler, respectively. 


The total moment of inertia and the mass of the each 
link with its counterweight are given by 
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II 

1 — 1 

M 

-r * 

If + 

hi’ 

[(X, - 

h)2 
2 ^ 

+ yp^J + 






[(^1 

0 

0 

1 

s a^) 

^ + (yp-^ sin a )2] 



I3 = 

+ 

1 

nw 

3 

[(X3 - 

i 3)2 
2 ^ 

+ y^^] + 


(D- 

-15) 


+ mj. 

[(^ 

3 

cos a 

3)^ + iyy 

-Rb sin cxb)^] 



and. 



mi = 

mi 

mA 











(D- 

-16) 




m3 = 

m,! 

3 






I'he radial distance and the angular position of the 
centroid for the crank and the coupler, respectively, are 
given by 


^1 

= i(x 3 

+ yp^) 

9i 

= 

(yp/^p) 


CM 

II 

+ y3^) 

93 

. -1 
= Sin 


the 

coupler extension 


(D- 17 ) 


about its centroid are given by 


mp 

~ 

^P ^P ’’^P ^ 

mR 

= 

^R ^R ^R 

Ip 

= 

1 T 2 

12 mp -^n 

Ir 

= 

1 rn r 2 

12 R "^R 

where P , b, t are 

the 

material 

of the extension 

and 

the suffix 

and right ends of 

the 

coupler. 


(D- 18 ) 

(D-19) 

mdth and thickness 
stand for the left 



I'he position of the centroid of the coupler is given 


by 




2 L R- 


(D-20) 


y2 = • 


ihe total mass and moment of inertia are given by 


^2 = “^2 ^R 


(D-21) 


Ip - ^2 ■*' ^ ^R \ 

(l)-22) 

I'he radial distsnce and angular position of the centroid of 
the coupler are given by 


rp = Xp 


cp = 0 


(D-25) 


Eqs. (D-15) to (D-17) and (D-21) to (D-23) give the 
parameters of the equivalent dynamical system needed to 
make the dynamic force analysis of the mechanism with 
counterweights. It must be noted that the effect of the 
holes in links or the pins for the joints and any enlargement 
in the links at the joints have been neglected in the above 
equations 5 presented for determining equivalent dynaiiJical 
s ys t em 1 

When Sector type counterweights are used at the 
coupler ends the eqs. (D-18) to (D.20) and eq. (D-22) are 


modified as follows 



(D-24) 


“l = ’' % -^L 

K 




-I 




(D-25) 


I 7-. — 2 Rq 


■R 


R R 


where 6 and 9 are the sector ajigles in radians of the 
left and right end counterweights, respectively- (Refer 
Pig. 5-11 in Chapter V.) 


, Ig . 4 sin ( ,/2) 

^2 [“^2 2 3 


4 r/2^, .. 

■L^ - - 3 r 


R. ) + mJl^ + 


/[m^ + 2^ + nij^] 


(D-26) 


, . sin ./2 2 

and 1 2 = 1 2 + Ijj '*■ (^2 '^ 3 — ^ ^ 


R 


4 


+ [X2 - (1^ + ^ 


2 2 


sin /2 p 

^^-R)]^ (I)-27) 

R 


where R^^ and R^ are the sector radii of the left and right 
end counterweights, respectively. 



APPENDIX B 


BUIiMARY OP COI-iPUi'ER PROGRANB 

A number of comparer programs in puRiRAlJ Iv’ have been 
written which are used to soItt-g the numerical examples in 
Chapter V. Piost of the programs are general and can be used 
in applications similar to the ones discussed in the present 
study for any planar mechanism. Ihese programs are described 
below. 

Program I ( S 1' Rnl N ) 

I'iiis program determines the deflection and stresses 
in any planar mechanism having one or more number of degrees 
of freedom. It also makes dynaniic force analysis. Ihe 
input data required is the link lengths and cross-sections, 
input crank speed of the mechanism and material properties 
of the links, i'his program consists of the following 
subroutines . 

(i) EKEM ; It determines the mass matrix and stiffness 
matrix of the element (link) in element oriented coordi- 
nates . 

(ii) iCLNANA : I'his program determines the angular position, 
absolute angular velocity and absolute angular accelera- 
tion of the various linms for a given constant input 


crank sneed . 



(iii) ELOjUj s It finds the element inertia load vector 
due to rigid body inertia of the links a.nd rlieir 
counterweights (if any) and the inertia load vector due 
to KED inertia oi the count erx-reights only in element 
oriented coordinates. 

(iv) yMPlV ; I'his program transforms the element oriented 
element matrices into system oriented element matrices 
and then assembles the element matrices into system 
matri ces . 

(v) EIG-SOL; I'his program solve the eigenvalue problem 
and can supply the eigenvalues and modal matrix for 
the case when the rigid body deg'ree of freedom is 
eliminated according to the procedure given in 
Chapter II and also when the crank input end is assumed 
instantaneously fixed as cantilever beam. I'his program 
requires a subroutine for the determination of eigen- 
A^-alues and eigenvectors of the dynamical matrix which 
is discussed in Appendix C. 

(vi) DBEAN; I'his program finds the displacement, velocity 
and acceleration vectors of the system at the end of 
a ceirtain chosen inter\ral by solving the equations 

of motion when initial displacement and velocity 
vectors ate given. 

(aO PROMAI; It is required to multiply matrices frequently 
during the aboA^-e process, i'his program is used for 
multiplication of matrices. 




Program II (OPTIMM) 

The program minimizes any iunction using PPP method 
of minimization. The input data is the initial design 
vector. It consists of the following subroutines. 

(i) This program is used to determine gradients of 
any function. In the present work it determines the 
gradient of the minimizing function, the maximum link 
stresses and the maximum shaking force. The former 
gradient is used in optimization process, whereas the 
latter ones are used to predict the vakues of latter 
mentioned functions during one dimensional minimization. 

(ii) PUN; This is a program in which the objective function 
and the various constraints are determined with the help 
of other subroutines. It also formulates the penalty 
function which is to be minimized. 

(iii) MIKGfOI; This program finds the optimum step length in 
a certain move direction. It first brackets the 
minimum within a certain interval and then the interval 
is reduced using Golden-section method to obtain the 
optimum step length. 

Profiram III (DIN' AN) 

This program makes the dynamic force analysis using 
the explicit eCLUations given in Appendix D. The input data 
is the length, mass, moment of inertia, polar coordinates of 
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of the centroid of each link. I'he input crank speed and 
the density of the material of the links. 

Pr ogram 17 (klNSYl) 

It is used to S 5 mthesi 2 '.e a planar four bar mechanism 
for rigid bodj guidance. I'he inpu'c data is the two angular 
relative positions of xhe coupler and two that of the 
follower. I'he rigid body positions and the corresponding 
crank rotations are already furnished to the program. 

Program V (ELPESI) 

Piiis program determines xhe mass, moment of inertia, 
position of the centroid of the links when counterweights 
are attached to them. I'he input data is the length, v?idth 
and thickness of each link and the radii of counterweights 
and. their angular positions, 'i'he density of the material 
is also supplied. 
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